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ABSTRACT 


The  Fast  Fourier  Transform  (FFT)  has  dramatically  changed 
the  signal  processing  field  by  allowing  the  transformation 
of  data  from  the  time  domain  to  the  frequency  domain  at 
speeds  rivaling  and  in  many  cases  exceeding  real  time.   The 
output  spectrum  obtained  from  the  FFT  depends  on  the  ampli- 
tudes and  the  phase  coherence  of  the  sinusoids  that  comprise 
the  input  waveform.   It  is  the  purpose  of  this  investigation 
to  modify  the  FFT  so  that  it  will  rely  more  heavily  on  the 
phase  coherence  of  these  sinusoids  rather  than  on  their  am- 
plitude.  In  this  way  a  second  indicator  of  signal  presence 
or  absence  would  be  available  to  add  more  reliability  to  the 
estimation  process.   The  study  proposes  a  method  for  reduc- 
ing the  dependence  of  the  output  coefficients  on  input  ampli- 
tudes and  investigates  its  merits  with  several  types  of  data. 
A  method  of  displaying  the  phase  coherence  indicators  is 
presented. 
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I.   INTRODUCTION 

Signal  processing  deals  with  the  problem  of  detecting 
signals  in  the  presence  of  noise.   If  the  noise  is  station- 
ary Gaussian  with  a  mean  value  of  0,  and  if  the  signal  is 
absolutely  stable,  then  it  is  theoretically  possible  to  dif- 
ferentiate between  the  signal  and  the  noise,  no  matter  how 
small  the  signal  to  noise  ratio.   By  integrating  long  enough 
(coherent  integration)  or  averaging  enough  transforms  (inco- 
herent integration)  the  noise  will  eventually  average  out  to 
zero  or  achieve  a  state  of  zero  variance.   This  is  the  abso- 
lute ideal  and  is  never  found  in  practice.   The  noise  is 
rarely,  if  ever,  purely  Gaussian,  nor  is  it  absolutely  sta- 
tionary.  Signals  are  not  absolutely  stable,  nor  do  they 
exist  for  infinite  periods  of  time,  thereby  disallowing  in- 
finite integration  or  averaging  times.   Because  of  these 
real-world  constraints  there  are  limitations  to  the  effec- 
tiveness of  all  signal  processing  methods  when  applied  to 
practical  problems. 

The  Fast  Fourier  Transform  (FFT)  is  one  of  the  most  com- 
mon and  effective  methods  by  which  data  can  be  transformed 
from  the  time  domain  to  the  frequency  domain.   This  trans- 
formation produces  a  frequency  spectrum  that  gives  the  rela- 
tive power  in  the  various  frequencies  that  comprise  the  time 
domain  data.   Due  to  the  practical  limitations  of  signal 
processing  stated  above  this  method  of  signal  processing  is 
limited  in  effectiveness.   It  is  the  purpose  of  this  study 
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to  try  to  develop  another  transformation,  similar  to  that  of 
the  FFT,  that  will  produce  a  frequency  spectrum  that  is  more 
dependent  on  the  phase  coherence  of  the  different  frequency 
components  that  constitute  the  time  domain  data  and  less  de- 
pendent on  their  power  content.   The  spectrum  analyst  will 
then  have  a  second  indicator  to  aid  him  in  making  his  deter- 
mination of  signal  presence  or  absence.   Hopefully  this  will 
result  in  more  reliable  estimates  and  allow  signal  detection 
in  some  situations  that  have  previously  been  beyond  the  phys- 
ical capabilities  of  the  FFT. 

Part  1  of  this  investigation  developes  a  procedure  for 
obtaining  a  numerical  indicator  of  phase  coherence  from  the 
FFT.   The  procedure  involves  weighting  the  amplitudes  of  the 
components  that  add  to  produce  each  of  the  FFT  output  coeffi- 
cients equally  so  the  resulting  spectrum  will  rely  more 
heavily  on  the  phase  relationships  between  additive  components 
than  on  their  amplitudes.   Part  2  compares  the  merits  of  the 
procedure's  signal  enhancement  capabilities  with  that  of  the 
FFT.   Part  3  summarizes  these  comparisons  and  Part  4  draws 
some  conclusions  about  the  differences  between  the  Modified 
FFT's  spectrum  and  that  of  the  FFT. 
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II.   PART  1  DEVELOPMENT 

A.   GENERAL 

The  Fast  Fourier  Transform  (FFT)  is  an  algorithm  for  com- 
puting the  Fourier  Coefficients  of  a  discretely  sampled  wave- 
form X„(t)  where  the  subscript  T  implies  X(t)  to  be  finite  in 
length  (of  length  T)  and  also  periodic  in  t  with  period  T. 
fpEl/T  is  the  fundamental  frequence  of  the  FFT. 

Appendix  A  gives  a  brief  description  of  the  Discrete  and 
Fast  Fourier  Transforms.   The  flow  graph  of  a  16  point  FFT 
along  with  some  of  the  notation  to  be  used  throughout  this 
paper  is  presented  in  Figure  1-1.   The  term  "level"  will  be 
used  to  denote  the  resultant  vector  formed  after  each  pass 
(matrix  multiplication)  of  the  algorithm.   Level  (0)  refer- 
ences the  sampled  version  of  the  input  waveform  in  its  natu- 
ral order.   Level  (1)  references  the  input  vector  after  bit 
reversal  has  taken  place.   Level  (2)  is  the  intermediate  re- 
sult (or  vector)  resulting  from  the  first  pass  (matrix  multi- 
plication) through  the  FFT  algorithm.   The  number  2  is  used 
to  imply  that  two  elements  of  the  input  vector  are  involved 
in  the  computation  of  each  element  of  the  intermediate  vector 
at  level  (2).   The  notation  is  similar  for  the  remainder  of 
the  algorithm.   A  circle  with  a  +  or  -  in  it  implies  addition 
or  subtraction  is  to  take  place.   A  circle  with  a  1,2,3,... 
in  it  implies  multiplication  by  W^ ,  W^ ,  W\  .  . .  where 
W  E  exp(-j2Tr/N). 
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The  magnitudes  of  the  output  coefficients  from  the  FFT 
depend  on  the  amplitude  of  the  input  signal  and  its  phase  co- 
herence throughout  the  transformation.   When  the  FFT  algo- 
rithm is  viewed  as  a  sequence  of  vector  (or  phasor)  rotations 
and  additions,  it  is  the  phase  coherence  of  a  signal  that 
determines  how  the  intermediate  vectors  will  add  in  produc- 
ing each  output  coefficient  of  the  FFT. 

Each  pass  through  the  algorithm  consists  of  a  sequence  of 
vector  rotations,  followed  by  a  vector  addition.   As  previous- 
ly stated,  the  length  or  magnitude  of  each  output  coefficient 
is  a  function  of  the  magnitude  of  the  intermediate  vectors 
and  of  the  manner  in  which  they  add  to  make  up  each  output 
coefficient.   The  term  phase  coherence  is  used  to  describe 
the  net  effect  of  the  phase  relationships  throughout  the  en- 
tire transform.   A  pure  sinusoid  would  have  phase  coherence 
throughout  the  entire  transform,  whereas  the  phase  relation- 
ship for  noise  would  be  completely  random  and  incoherent . 
It  is  phase  coherence  that  allows  the  intermediate  vectors 
of  a  pure  sinusoid  to  add  end-to-end  in  a  straight  line, 
whereas  the  intermediate  vectors  for  noise  add  at  random 
phase  angles,  thereby  producing  output  coefficients  with  an 
amplitude  less  than  the  maximum  value  they  could  achieve  if 
they  were  to  add  in  a  straight  line. 

To  give  a  better  understanding  of  the  concept  of  phase 
coherence  consider  the  vector  additions  in  Figure  1-2.   The 
resultant  vector  may  be  viewed  as  the  output  vector  of  the 
FFT  and  its  component  vectors  as  the  intermediate  vectors  of 
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Vector  addition  analogy  of  FFT 

(a)  Incoherent  vector  additions 

(b)  Coherent  vector  additions 

(c)  Incoherent  vector  additions,  same  as  (a) 
but  with  intermediate  components  normalized 
to  same  length. 


Figure  1 -2 
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the  FFT  that  have  been  rotated  and  added  to  form  the  result- 
ant or  output  vector.   Figure  l-2a  could  be  viewed  as  inco- 
herent noise  with  no  particular  order  in  the  manner  in  which 
the  component  vectors  add  to  produce  the  resultant  vector. 
Conversely,  Figure  l-2b  shows  a  resultant  vector  of  what  may 
be  interpreted  as  the  result  of  the  transformation  of  a  pure 
sinusoid  with  zero  phase  or  frequency  shift  throughout  the 
time  record  being  transformed.   Each  component  adds  to  the 
others  in  a  straight  line  and  at  the  phase  angle  of  the  input 
sinusoid.   The  vector  in  2a  achieves  its  magnitude  not 
through  phase  coherence  but  simply  because  one  of  its  com- 
ponent vectors  dominates  the  others  in  magnitude.   If  each 
of  its  component  vectors  were  restricted  to  having  unit 
length  while  still  maintaining  its  same  phase  angle,  then 
the  magnitude  of  the  resultant  vector  would  be  forced  to  de- 
pend on  phase  coherence  (the  ability  of  the  component  vectors 
to  add  in  a  straight  line)  throughout  the  transform  and  not 
on  the  magnitude  of  any  of  its  component  vectors.   This  sit- 
uation is  depicted  in  Figure  l-2c. 

When  the  output  of  the  discrete  correlation  performed  by 
the  FFT  is  squared,  the  resulting  spectrum  is  proportional 
to  the^ power  contained  in  the  original  waveform  at  each  of 
the  correlating  frequencies  n/T  (n=0 , 1 , 2 , . . . , N/2-1 )  for  real 
valued  input  waveforms  and  is  referred  to  as  the  power  spec- 
tral density  of  the  input  waveform.   The  association  of  the 
output  coefficients  of  the  FFT  is  strongly  related  to,  and 
dependent  on,  the  power  of  the  input  waveform  at  its  various 
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frequencies,  or  more  explicitly,  at  the  correlation  frequen- 
cies of  the  FFT  algorithm. 

It  would  seem  plausible  that  if  the  magnitude  of  the  in- 
termediate vectors  could  be  removed  then  the  output  spectrum 
would  rely  more  on  the  sinusoidal  nature  of  the  input  wave- 
form than  on  the  power  of  the  input  sinusoid.   Or  conversely, 
the  ability  of  large  vectors,  due  to  incoherent  noise,  to 
corrupt  or  mask  a  coherent  sinusoid  would  be  reduced.   The 
objective  is  to  develop  an  indicator  of  the  presence  of  a 
signal  through  phase  coherence  with  as  little  emphasis  as 
possible  on  signal  or  noise  power,  an  indicator  which  could 
be  used  in  conjunction  with  the  FFT  to  aid  in  estimating  sig- 
nal presence  or  absence. 

B.   PREMISE 

The  basic  premise  of  this  study  is  to  modify  the  FFT 
algorithm  in  such  a  manner  that  the  output  coefficients  give 
"  an  indication  of  the  relative  phase  coherence  of  the  sinu- 
soids that  make  up  the  input  waveform.   To  test  the  modified 
transform  to  see  if  it  meets  the  premise,  sinusoids  that 
complete  an  integral  number  of  cycles  across  the  time  record 
will  be  considered,  that  is,  sinusoids  that  are  integral 
multiples  of  the  fundamental  frequency  of  the  transform  fo. 
These  are  of  course  the  correlating  sinusoids  of  the  FFT  algo- 
rithm.  It  will  be  required  to  show  that  the  phase  angle  of 
the  output  coefficient  corresponding  to  a  pure  sinusoidal 
input  at  any  frequency  nf o (n=l , 2 , . . . , N/2-1)  reflects  the 
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phase  angle  of  the  input  sinusoid.   This  must  be  true  if  the 
intermediate  vectors  are  to  add  in  a  straight  line  for  a 
phase  coherent  sinusoid. 

If  the  resulting  output  coefficient  representing  an  in- 
put sinusoid  of  frequency  nfo  reflects  the  phase  angle  of  an 
input  sinusoid  at  these  frequencies,  then  the  premise  has 
been  satisfied.   Sinusoids  at  frequencies  other  than  integral 
multiples  of  fo  are  seen  by  the  correlation  process  as  a  sum- 
mation of  sinusoids  at  the  frequencies  nfo.   That  is,  the 
FFT  output  does  not  reflect  the  phase  angle  of  the  input 
sinusoids  at  frequencies  other  than  integral  multiples  of  fo 
but  instead  approximates  these  sinusoids  by  sinusoids  that 
are  at  integral  multiples  of  fo.   Therefore  the  phase  angles 
obtained  do  not  reflect  the  phase  angle  of  the  input  sinu- 
soid, but  instead  reflect  the  phase  angles  of  the  correlating 
sinusoids  that,  when  summed,  will  approximate  the  input  wave- 
form.  Therefore  it  is  meaningful  to  test  the  premise  on 
input  sinusoids  that  are  integral  multiples  of  fo. 

C.   IMPLEMENTATION  OF  THE  MODIFICATION 

The  output  coefficients  of  the  FFT  are  in  general  com- 

X 

plex  numbers  of  the  form  c=a+jb  with  magnitude  |c|  =  (a^+b^)^ 
and  phase  angle  0=Tan~ ^ (b/a) .   When  an  output  coefficient  of 
the  FFT  is  viev/ed  with  the  objective  of  determining  the  phase 
angle  of  its  corresponding  input  sinusoid,  the  phase  angle 
obtained  from  0=Tan~'(b/a)  is  the  average  difference  in  phase 
angle  between  the  input  sinusoid  and  a  reference  cosine.   If 
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a  sine  wave  is  used  for  the  reference  sinusoid  then  the  phase 
angle  e=Tan~^(a/b)  must  be  used.   In  addition,  one  must  know 
whether  the  correlating  sinusoids  were  of  the  form  cos(2Trnk/N) 
-jsin(2iTnk/N)  or  cos(2Tynk/N)+jsin(27rnk/N) .   Either  form  re- 
sults in  the  same  magnitude  but  with  a  different  phase  angle. 
When  6  is  computed  as  Tan~Ht)/a),  changing  the  sign  of  the 
sine  wave  correlation  results  in  a  reflection  of  the  phase 
angle  across  the  real  axis.   When  9  is  computed  as  Tan~''(a/b), 
the  sign  change  results  in  a  reflection  across  the  imaginary- 
axis.   The  phase  angle  has  physical  meaning  if  one  recognizes 
what  the  reference  sinusoid  was  and  whether  exp(-j2Trnk/N)  or 
exp( j2TTnk/N)  was  used  for  the  correlation  process. 

After  each  pass  of  the  FFT  algorithm  (at  each  level),  the 
intermediate  coefficients  can  likewise  be  represented  in  the 
same  manner  as  the  output  coefficients.   These  intermediate 
coefficients  are  analogous  to  the  intermediate  vectors  re- 
ferred to  in  Figure  1-2.   The  object  will  be  to  make  these 
intermediate  coefficients  unit  vectors,  vectors  of  magnitude 
1.   In  this  way  the  phase  relationships  of  these  vectors  will 
dominate  the  length  of  any  given  output  coefficient  with  as 
little  emphasis  on  the  magnitude  of  these  vectors  as  possible. 
The  method  utilized  will  be  simply  to  divide  each  coeffi- 
cient by  its  magnitude.   This  process  will  be  referred  to  as 
normalization  since  it  normalizes  the  magnitude  of  the  coef- 
ficients to  a  value  of  1  and  in  no  way  alters  the  ratio  of 
the  real  to  imaginary  part  of  each  coefficient,  that  is,  the 
phase  angle  of  each  coefficient  is  retained. 
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The  problem  then  becomes  one  of  determining  where  the 
normalization  should  take  place  in  the  FFT  algorithm.   The 
following  properties  of  the  FFT  butterfly  will  be  used  to 
make  this  determination  for  real  valued  input  waveforms: 

(1)  All  of  the  coefficients  at  the  2  level  are  real 
valued  since  each  coefficient  at  this  level  is  either  the 
sum  or  difference  of  2  of  the  original  real  valued  input  co- 
efficients. 

(2)  Half  of  the  coefficients  at  the  4  level  are  complex 
niimbers,  namely  the  odd  numbered  coefficients  XJ^  (k=l,  3,  5,  .  .  .  , 
N-l)  and  these  can  have  any  phase  angle,  the  rest  are  real 
valued. 

(3)  The  odd  numbered  coefficients  at  the  4  level  affect 
only  the  odd  numbered  output  coefficients  (those  associated 
with  fo  and  its  odd  harmonics),  and  similarly,  the  even  num- 
bered coefficients  affect  only  the  even  numbered  output  co- 
efficients (those  associated  with  the  even  harmonics  of  fo)- 

(4)  Above  the  4  level  an  increasing  number  of  coeffi- 
cients become  complex,  for  example  at  the  8  level  the  complex 
coefficients  are:  X®  (k=l , 2 , 3 , 5 , 6 , 7 , 9 , 10, 11 , . . . , N-l ) .   At  the 
16  level  the  complex  coefficients  are:  X^^  (k=l , 2 , 3 , 4 , 5 , 6 , 7 , 
9,10,11,12,13,14,15, .. . ,N-1),  etc. 

If  the  input  coefficients  are  normalized,  the  only  pos- 
sible phase  angles  that  they  can  have  are  0  or  180  degrees 

(each  with  a  magnitude  of  1)  since  they  are  all  real  numbers. 

N 
With  N  input  coefficients  then  there  will  only  be  2   possible 

normalized  input  vectors  (if  the  value  zero  is  allowed  then 
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N 
there  will  be  3  possible  input  vectors).   This  results  in  a 

finite  number  of  output  phase  angles  and  thus  by  contradic- 
tion prohibits  the  vectors  from  adding  in  such  a  manner  as 
to  reflect  all  possible  input  phase  angles. 

The  same  argument  is  true  for  the  2  level  where  the  co- 
efficients again  are  real  numbers.   In  fact  the  number  of 
possible  output  phase  angles  is  less  than  when  the  coeffi- 
cients were  normalized  at  the  input.   The  number  of  differ- 

N 
ent  vectors  is  still  2   at  the  2  level  but  there  is  one  less 

pass  between  the  2  level  and  the  output  than  between  the  in- 
put and  the  output,  thereby  reducing  the  number  of  passes  in 
which  coefficients  are  combined. 

The  problem  of  determining  where  to  implement  this  nor- 
malization procedure  is  greatly  simplified  if  one  remembers 
that  the  purpose  of  this  modification  is  to  achieve  vector 
additions  of  the  intermediate  coefficients  so  that  these  vec- 
tors will  add  in  such  a  way  that  the  resultant  vector  is 
oriented  in  the  direction  of  the  input  waveform's  phase.   In 
general  that  phase  will  be  measured  with  respect  to  a  cosine 
waveform  throughout  this  study. 

What  the  above  suggests  is  that  a  coefficient  should  be 
complex  before  it  is  normalized  in  order  not  to  restrict  the 
•total  number  of  phase  angles  achievable  by  the  output  coeffi- 
cients . 

The  4  level  is  the  first  level  where  complex  numbers 
arise.   As  previously  stated  half  of  the  coefficients  affect 
only  the  output  coefficients  related  to  even  harmonics  of  fo 
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At  this  level  the  transform  can  be  split  up  into  two  distinct 
and  independent  parts,  one  for  the  odd  harmonics  and  one  for 
the  even  harmonics.   The  real  coefficients  are  the  building 
blocks  of  the  even  harmonics,  by  the  previous  argument  of  a 
finite  number  of  output  phase  angles  it  will  not  be  possible 
to  normalize  these  coefficients  and  obtain  all  possible 
phase  angles  at  the  output. 

The  complex  coefficients  at  this  level  are  the  interme- 
diate coefficients  that  make  up  the  output  coefficients  re- 
lating to  the  odd  harmonics  of  fo.  If  these  coefficients 
are  normalized  the  output  coefficients  corresponding  to  the 
odd  harmonics  of  fo  will  reflect  the  phase  angle  of  any  in- 
put sinusoid  at  the  frequencies  f o , 3f o , 5f o , . . . , (N/2-l)f o . 
Demonstration: 

Consider  the  bit  reversed  form  of  a  16  point  transform 
of  some  function  x(t)  and  the  FFT  butterfly  for  these  coef- 
ficients up  to  the  4  level  in  Figure  1-3.   As  the  notation 
indicates  each  coefficient  at  the  4  level  is  made  up  of  four 
elements  of  the  input  vector.   Each  grouping  of  four  coeffi- 
cients at  the  4  level  can  be  rewritten  as: 
X'*(4j)    =  x(0+j)  +  x(N/2+j)  +  x(N/4+j)  +  x(3N/4+j) 
X'*(4j+1)  =  x(0+j)  -  x(N/2+j)  +  j[x(N/4+j)  -  x(3N/4+j)] 
X'*(4j+2)  =  x(0+j)  +  x(N/2+j)  -  [x(N/4+j)  +  x(3N/4+j)] 
X'*(4j  +  3)  =  x(0+j)  -  x(N/2+j)  -  j[x(N/4+j)  -  x(3N/4+j)] 
where  j=0, 1 , 2 , . . . ,N/4-l .   This  is  a  valid  generalization  for 
any  N.   The  coefficients  in  question  are  the  complex  coeffi- 
cients and  only  these  coefficients  will  be  considered.   Since 
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(0) 
x(0) 
x(1) 
x(2) 
x(3) 
x(l|) 
x(5) 
x(6) 

x(7) 
x(8) 

x(9) 
x(10) 
x(11) 
x(12) 

x(13) 
x(1l4.) 
x(l5) 


(1) 
x(O),^— — ^x(0 
x(8)^-"^^x(0 

x(2)^.;— -^x(2 
x(10)-^^~^x(2 

X(6)^;;— — r  x(6 

x(l)^.— — =.x(1 
x(9)-^^  x(1 
x(^)^— -^x(5 
x(13)-^^x(5 
x(3)^^ ^x(3 


x(11) 
x(7) 


x(l5) 


x(3 
x(7 


x(7 


+x(8)+x(l4.)+x(12) 

-x(8)  +  j(x(i|)-x(12)) 

+x(8)-(x(l4)+x(12)) 

-x(8)-j(x(l|)-x(12)) 

+x(10)+x(6)+x(1l4.) 

-.x(10)  +  j(x(6)-x(1i;)) 

+x(10)-(x(6)+x(1I+)) 

-x(10)-j(x(6)-x(1l4.)) 

+x(9)+x(5)+x(13) 

-x(9)+j(x(5)-x(13)) 

+x(9)-(x(5)+x(13)) 

-x(9)-j(x(5)-x(13)) 

+x(11)+x(7)+x(l5) 

-x(11)+j(x(7)-x(15)) 

+x(11)-(x(7)+x(l5)) 

-x(11)-j(x(7)-x(l5)) 


The  FFT  algoritm  for  a  16  point  transform  up  to  the  i^.  level 


Figure  1-3 
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normalization  is  simply  the  division  of  any  coefficient  by 
its  magnitude,  then  the  normalization  factor  (the  magnitude 
of  the  complex  coefficients)  is: 

{[x(j)  -  x(N/2+j)]'  +  [x(N/4+j)  -  x(3N/4+j)]M* 
The  input  sinusoid  can  be  generalized  as: 
x(t)  =  Asin(2TTnfot  +  6)  where  n=l,  3,  5,  .  .  .  ,N/2-l,  fo=l/T 
and  6  is  arbitrary  -•iT_<6£Tr. 

The  sampled  version  of  x(t),  namely  x(k)  (k=0, 1, 2 , . . . ,N-1) 
is:   x(k)=Asin(2TTnf okAT  +  6)  where  AT  is  the  time  increment 
between  samples. 
Now:   T=NAT  and  f o=l/T=l/ (NAT) ,  then: 

x(k)=Asin(2TTnk/N  +9). 
Therefore  the  divisor  of  the  coefficients  can  be  written, 

{[Asin(2^nj/N+6)-Asin(2TTn(N/2+j)/N+e)]2 

+  [Asin(27Tn(N/4+j)/N+e)-Asin(2Trn(3N/4+j)/N+e)]M^ 
letting  C=2Trnj/N  and  B=2TTn,  then  the  divisor  can  be  rewritten: 
A{ [sin(C+e)-sin(C+B/2+e)]2  +  [sin(C+B/4+e )-sin(C+3B/4+e ) ] M^ 
Now,  sin(x+y )=sin(x)cos(y )+sin(y )cos(x) 
Expanding  the  above, 

A{[sin(C)cos(e)+sin(e)cos(C)-sin(C)cos(B/2+e)-sin(B/2+e)cos(C)]^ 
+[sin(C)cos(B/4+e)+sin(B/4+e)cos(C) 

X 

-sin(C)cos(3B/4+e)-sin(3B/4+6)cos(C)]M^ 
Now ,  cos ( x+y )=cos(x)cos(y)-sin(x)sin(y) 
Expanding  once  more. 
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A{ [sin(C)cos(e)+sin(0)cos(C)-sin(C)cos(B/2)cos(e) 

+sin(C)sin(B/2)sin(e)-cos(C)sin(B/2)cos(0) 

-cos(C)sin(e)cos(B/2)]2    +    [sin(C)cos(B/4)cos( 9) 

-sin(C)sin(B/4)sin(0)+cos(C)sin(B/4)cos(e) 

+cos(C)sin(e)cos(B/4)-sin(C)cos(3B/4)cos(e) 

+sin(C)sin(3B/4)sin(9)-cos(C)sin(3B/4)cos(e) 

-cos(C)sin(e)cos(3B/4)]2}* 

Now  B= 2 Tin    (n=l,3,5,  .  .  .  ,N/2-l) 

=  2Tr,6TT,  lOlT,  14lT,  .  .  . 

then   B/2=it,3tt,  5tt,7tt,  .  .  . 

and   sin(B/2)=0 

cos(B/2)=-1 

B/4=tt/2,3tt/2,5tt/2,7tt/2,  .  .  . 
sin(B/4)  =  ±l=(-l)f<"-^l)/2]-M 

cos(B/4)=0 

3B/4=3TT/2,9TT/2,15Tr/2,2lTT/2,  .  .  . 

sin(3B/4 )  =+!=(-!)  ^'^"^^^/^ 

cos(3B/4)-0 
Making  the  above  substitutions: 

A{ [sin(C)cos(e)+sin(e)cos(C)+sin(C)cos(e)+cos(C)sin(e) ] ^ 
+[+sin(C)sin(e)±cos(C)cos(e)+sin(C)sin(G)±cos(C)cos(0) ] M^ 
=A{ [2sin(C)cos(e)+2sin(e)cos(C)] ^+ [+2sin(C)sin( G ) 
±2cos(C)cos(e)]M*   =    A{4sin2(C)cosHe)  +  8sin(C)cos(e)sin(6)cos(C) 
+4sin^(e)cos^(C)+4sin^(C)sin^(e)-8sin(C)sin(e)cos(C)cos(e) 
+4cos^(C)cos^(e)}*   =    2A{sin^(C) [cos^(e)+sin^(e)] 
+cos(C)  [sin^(e)+cos^(0)]  }^'   =    2A{sin^  (C)+cos^  (C) }  ="    =    2A 
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Or,  the  divisor  at  the  4  level  for  the  complex  coefficients 
is  dependent  only  on  the  amplitude  of  the  input  sinusoid;  it 
is  independent  of  N,n,j,  or  9,   This  is  equivalent  to  divid- 
ing all  of  the  coefficients  by  a  constant  for  any  given  in- 
put sinusoid  of  the  form  x(t  )=Asin(2TTnf  ot  +  O)  and  in  no  way 
alters  the  phase  angle  of  the  output  coefficients.   Its  net 
effect  is  that  of  scaling  the  input  sinusoid  such  that  it  is 
replaced  by: 

x(t  )  =  0.  5sin(2iTnf  ot  +  9)   for  any  amplitude  A. 

It  can  likewise  be  shown  that  a  similar  normalization  at 
the  8  level  results  in  division  of  the  coefficients  by  4A, 
at  the  16  level  the  divisor  is  8A  and  in  general  the  divisor 
will  be  (L/2)A  where  L  is  the  level  at  which  the  normaliza- 
tion was  performed. 

The  second  half  of  the  transformation  can  not  be  normal- 
ized at  the  4  level  for  reasons  already  stated.   At  the  8 
level  the  coefficients  once  more  divide,  that  is,  the  trans- 
formation can  be  broken  down  into  4  independent  transforma- 
tions, each  part  being  related  to  one-fourth  of  the  output 
coefficients.   The  intermediate  coefficients  1,5,9,13,.... 
make  up  their  corresponding  output  coefficients,  likewise 
3,7,11,15,...  make  up  their  corresponding  output  coefficients, 
These  coefficients  offer  no  problem  since  they  were  already 
complex  at  the  4  level.   Coefficients  0,4,8,12,...  are  still 
real  valued  at  the  8  level  and  do  not  lend  themselves  to 
normalization  at  this  point.   Coefficients  2,6,10,14,... 
were  not  complex  at  the  4  level  but  are  complex  at  the  8 
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level.   A  similar  proof  to  the  one  used  for  normalization  of 
the  intermediate  coefficients  of  the  odd  harmonics  of  fo  at 
the  4  level  can  be  used  to  show  that  normalizing  the  coeffi- 
cients 2,6,10,14,...  at  the  8  level  does  not  alter  the  phase 
angle  of  the  output  coefficients  corresponding  to  these  fre- 
quencies for  sinusoidal  inputs  that  are  integral  multiples 
of  fo.   Similarly  at  the  16  level  the  intermediate  coeffi- 
cients that  make  up  the  frequencies  nfo  (n=4, 12, 20, 28, . . . ) 
become  complex  and  can  therefore  be  normalized  without  alter- 
ing the  phase  angle  of  the  output  coefficients  corresponding 
to  these  frequencies. 

This  process  continues  until  one  reaches  the  output  of 
the  FFT  at  which  time  all  of  the  coefficients  (except  F(0), 
the  D.C.  term)  have  become  complex  and  therefore  normaliza- 
tion could  have  been  performed.   The  idea  of  normalizing  dif- 
ferent coefficients  at  different  levels  produces  one  serious 
drawback.   Normalization  fixes  the  maximum  magnitude  an  out- 
put coefficient  can  achieve.   This  is  due  to  the  fact  that 
every  intermediate  coefficient  somewhere  in  the  FFT  algorithm 
has  been  replaced  with  a  magnitude  of  1.   Knowing  where  the 
normalization  occurred  one  then  knows  how  many  more  passes 
will  be  performed  before  the  algorithm  is  complete.   Since 
each  pass  consists  of  rotating  and  then  adding  coefficients, 
the  number  of  vector  additions  remaining  per  output  coeffi- 
cient is  then  equal  to  the  number  of  remaining  passes.   Since 

the  coefficients  have  magnitude  1,  they  can  achieve  a  magni- 

p 
tude  at  the  output  no  larger  than  2  where  P  is  the  number  of 
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remaining  passes  to  be  performed.   For  instance,  in  a  32 
point  FFT,  if  the  intermediate  coefficients  that  make  up  the 
odd  harmonics  of  fo  are  normalized  at  the  4  level  there  will 
be  three  remaining  passes  yet  to  be  performed.   At  the  8 
level  their  maximum  magnitude  is  2,  at  the  16  level  it  is  4 
and  at  the  32  level  (the  output  in  this  case)  their  magnitude 
can  be  8.   The  only  way  the  value  8  can  be  achieved  is  if  the 
input  was  in  fact  a  pure  sinusoid  at  a  frequency  nfo  (n=l,3, 
5,...).   It  is  important  to  note  that  if  one  knows  how  many 
points  (N)  were  in  a  transformation,  and  what  frequency  a 
particular  output  coefficient  is  related  to,  that  is,  if  one 
knows  N,  then  the  maximum  value  that  coefficient  can  achieve 
is  also  known.   From  this  information  one  may  be  able  to 
make  a  prediction  as  to  the  probability  of  the  existence  of 
a  signal  at  the  frequency.   If  one  also  considers  the  magni- 
tude of  this  coefficient  from  more  than  one  transform,  then 
his  prediction  becomes  more  reliable.   This  in  itself  could 
be  a  powerful  tool  in  signal  processing  work  since  the  out- 
put coefficients  are  independent  of  the  amplitude  of  the  in- 
put signal. 

The  serious  drawback  occurs  when  one  analyzes  a  spectrum 
where  the  intermediate  coefficients  have  been  normalized  at 
different  levels.   It  is  a  simple  procedure  to  multiply  each 
output  coefficient  by  the  number  of  the  level  at  which  it 
was  normalized.   In  the  previous  example  a  32  point  transform 
was  considered  and  it  was  shown  that  the  maximum  magnitude 
that  an  output  coefficient,  dependent  on  normalization  of 
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its  intermediate  coefficients  at  the  4  level,  could  achieve 
was  8.   Then  8  times  4  =  32.   Likewise  an  output  coefficient 
dependent  on  normalization  of  its  intermediate  coefficients 
at  the  8  level  can  achieve  a  magnitude  of  4.   Then  4  times 
8  =  32.   In  general,  the  magnitude  of  the  output  coefficients 
can  be  made  a  constant  (namely  N)  for  an  input  sinusoid  of 
the  form  Asin(27rnf  ot  +  0)  (n=l,  2,  3,  .  .  .  ,N/2-l)  independent  of 
A,  n,  and  6.   And  the  phase  angle  9  can  be  obtained  from  the 
output  coefficients. 

When  the  input  waveform  consists  of  signal  plus  noise, 
the  multiplication  process  leads  to  erroneous  results.   To 
best  explain  this  problem,  consider  the  output  coefficient 
which  is  dependent  on  normalization  at  the  N    level.   This 
will  always  be  the  frequency  (N/4)fo  since  this  frequence  co- 
efficient does  not  become  complex  until  after  the  last  pass 
of  the  FFT  algorithm  has  been  performed.   This  coefficient 
will  always  have  a  magnitude  of  1  after  normalization  and  it 
will  have  a  magnitude  N  after  it  has  been  multiplied  by  its 
normalization  level,  which  is  always  N  in  this  case.   That 
is,  the  frequency  (N/4)fo  will  always  have  the  magnitude  N 
regardless  of  whether  there  was  signal  or  noise  at  that  fre- 
quency.  To  a  lesser  degree  this  same  problem  extends  back 
to  the  frequencies  that  were  normalized  at  the  N/2  level, 
specifically  (N/8)fo  and  (3N/8)fo.   In  general  the  frequen- 
cies that  are  normalized  late  in  the  transform  tend  to  domi- 
nate the  spectrum.   This  can  be  explained  if  one  considers 
that  the  multiplication  factor  was  simply  an  artificial  means 
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of  adding  amplitude  to  the  output  coefficients  so  that  all 
of  the  frequencies  could  achieve  the  same  maximum  value. 
Frequencies  that  have  their  components  normalized  early  in 
the  transform  rely  more  heavily  on  phase  coherence  to  achieve 
large  amplitudes.   The  frequencies  that  have  their  components 
normalized  late  in  the  transform  have  less  time  to  achieve 
large  amplitude  through  phase  coherence. 

One  solution  to  the  above  problem  is  to  compare  coeffi- 
cients that  are  dependent  upon  normalization  at  a  given  lev- 
el independent  of  the  remaining  coefficients.   This  has  the 
drawback  that  the  reliability  of  the  estimate  of  signal  pres- 
ence is  reduced  since  the  comparison  is  made  against  fewer 
coefficients.   In  addition  it  is  inconvenient  and  lacking  in 
continuity  since  at  best  every  other  coefficient  is  compared 
leaving  gaps  between  frequencies. 

A  second  solution  was  devised  in  which  all  of  the  coef- 
ficients that  depend  on  normalization  at  a  given  level  are 
summed.   Their  mean  value  is  determined  and  then  each  coef- 
ficient is  divided  by  this  mean  value.   This  has  the  effect 
of  giving  the  coefficients  dependent  on  each  level  of  normal- 
ization a  mean  value  of  1.   The  spectrum  can  then  be  recom- 
bined  with  the  effect  of  giving  the  entire  spectrum  a  mean 
value  of  1.   A  quick  example  will  illustrate  this  procedure. 

Consider  a  16  point  FFT  with  an  input  signal  consisting 
of  equal  amplitude  1  hertz  and  2  hertz  sinusoids  plus  noise. 
Figure  1-4  shows  the  input  sinusoids  and  the  sinusoids  plus 
noise.   Figure  1-5  shows  the  FFT  algorithm  and  Figure  1-6 

29 


Signal 


Signal  +  Noise 


0 
1 
2 

3 
k 
5 
6 

7 
8 
9 
10 
11 
12 
13 
1i^ 
15 


0.0 

1 .0898 
1.7071 
1.6310 

1.0 

0.2168 
•0.2929 
•O.32M1 

0.0 

0.32hk 
0.2929 
■0.2168 
•1.0 
■  1  .6310 
•1  .7071 
.1.0898 


0.1 
0.6 

2.0 

l.i^ 
1.2 

-0.3 

-1.0 

0.0 

-0.[|. 

o.S 
-o.k 
-0.5 

-0.9 

-1.9 

1.6 

-2.0 


1  hz  +  2  hz  sinusoids  +  noise 


Figure  1  -i| 


30 


X3 

Si 

5 

X! 

^- 

CM 

CA 

LA 

vO 

r- 

^ 

rA 

CA 

■LA 

rA 

o 

>JD 

rA 

O 

fA 

C^ 

tA 

nO 

XA 

O 

vO 

O^ 

C^ 

.d- 

,_ 

o 

J- 

rA 

rA 

CM 

^ 

-d- 

o 

• 

1 

• 

o 

• 

■LA 

« 

CM 

• 
rA 

• 

CA 

• 

CO 

O 

r- 

CA 

O 

CO 

CO 

o 

rA 

r-- 

O 

CO 

lA 

-d- 

o 

CM 

-=i- 

CO 

GO 

-d 

CM 

o 

,-t- 

XA 

c^ 

co 

-d 

_d 

CO 

c^ 

c— 

CO 

_d 

_d 

CO 

CT^ 

CT^ 

rA 

v£) 

CO 

rA 

o 

o 

rA 

CO 

sD 

CA 

O^ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•LA 

LA 

^— 

CM 

rA 

CM 

CM 

rA 

CM 

,— 

XA 

XA 

-^"^ 

■r-3 

'•-) 

•«-s 

•r-j 

•'-5 

••-5 

••-> 

'r-i 

•r-J 

•f-i 

•r-5 

•r-3 

vD 

+ 

+ 

1 

+ 

+ 

+ 

1 

1 

1 

+ 

1 

1 

<— 

C7^ 

CM 

CO 

CM 

CO 

rA 

5— 

CO 

CM 

CO 

CM 

C^ 

■ — ' 

-d- 

C^ 

C^ 

T- 

CM 

vO 

-X) 

CM 

cr^ 

r— 

.-t 

cr> 

r— 

rA 

rA 

CO 

r— 

r— 

CO 

rA 

rA 

c~ 

O 

o 

CA 

rA 

O 

CM 

O 

CO 

o 

-d- 

x£) 

CO 

O 

CM 

O 

rA 

CA 

o 

CM 

O 

-^ 

CM 

CM 

O 

lA 

-d 

XA 

O 

CM 

CM 

-d- 

O 

CM 

' 

1 

1 

1 

1 

' 

1 

1 

1 

UN 

u\ 

XA 

■LA 

r- 

r- 

r- 

r- 

CO 

K— 

^-~ 

CO 

r- 

f- 

r^ 

i>- 

•LA 

-d- 

-d- 

LA 

XA 

XA 

LA 

XA 

c^ 

O 

CM 

CM 

o 

O^ 

rA 

c^ 

c" 

T— 

O 

rA 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

^— 

^— 

Al 

CM 

^— 

,— 

-d- 

CM 

,— 

,— 

CM 

-d 

/--^ 

••-> 

•f-j 

•>-} 

•f-5 

"-i 

•"-J 

•<-3 

•r-3 

•r-j 

•"-5 

•1-5 

•r-5 

CO 

+ 

+ 

1 

+ 

t 

1 

+ 

+ 

+ 

1 

1 

1 

^ 

■LA 
•LA 
rA 

LA 
■LA 
fA 

LA 
■LA 
CA 

' 

XA 
XA 
CA 

rA 

CM 

O 

o 

CA 
CM 

OJ 

O 

vO 

O 

CM 

O 

v£i 

O 

CM 

O 

rA 

T- 

O 

r- 

rA 

O 

CM 

-d- 

o 

1 

rA 

CM 

CA 

1 

O 

1 

-d- 

CM 

O 

rA 

o 

o 

O 

rA 

O 

^— 

^_ 

^ 

vX) 

vO 

sO 

o 

O 

• 

• 

• 

• 

• 

• 

• 

• 

00 

CM 

CM 

CM 

5*" 

T— 

CM 

OJ 

^. s. 

•f-J 

•r-i 

••—i 

•<-} 

•1-3 

•""3 

••~i 

•<-5 

-d- 

+ 

1 

1 

+ 

+ 

1 

+ 

1 

O 

•LA 

xX) 

■LA 

CM 

-d- 

o 

^ 

'- 

■^ 

rA 

«r— 

C^ 

O 

0^ 

O 

o 

o 

1 

O 

CM 

CM 

^ 

CM 

7 

o 

rA 

o 

T 

s— 

CM 

^r- 

, , 

c<\ 

•LA 

CA 

^ 

vO 

-d 

vD 

vO 

%— 

^_ 

CM 

v£) 

o^ 

CJ- 

O 

o 

CM 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

* 

• 

• 

• 

• 

• 

' — 

o 

O 

o 

CM 

'" 

CM 

O 

AJ 
1 

^— 

o 

CM 

'^ 

o 

T- 

CM 
1 

CM 

^ 

r- 

-i 

CM 

O 

O 

^ 

o 

^X) 

vD 

XA 

r^ 

o 

-d 

XA 

O 

O 

^ 

o 

o 

^~ 

o 
1 

AJ 

O 

1 

r— 

^— 

O 

O 

O 
1 

^ — 

r- 

o 
1 

o 

C^J 

,.^ 

^ 

vD 

O 

-=f 

CM 

rA 

o 

o 

-d 

XA 

^ 

XA 

0^ 

0> 

vD 

O 

O 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

- — ' 

o 

O 

CM 

^— 

T~* 

O 

,— 

o 

o 

o 

O 

o 

o 

T— 

T- 

CM 

- 

1 

1 

1 

1 

1 

1 

1 

1 

rH 

I^H 

0) 
> 

0^ 

0) 

o 

o 

CM 

rA 

-d 

\r\ 

vO 

C-- 

CO 

o^ 

o 

T— 

CM 

rA 

-d 

XA 

»-l 

o 

s— 

*" 

'" 

^— 

T— 

<— 

id 

XA 

-p 

1 

•H 

T— 

U 

O 

Q) 

W) 

^ 

rH 

0 

< 

W 

•rH 

H 

fe 

t^ 

k. 

&H 
fe 
fo 


«) 

•H 
<;^ 
•H 
■C 
O 


• 

N 

N 

N 

N 

N 

N 

N 

• 

Si 

xi 

^ 

:& 

Si 

X! 

X! 

Q 

^ 

CM 

fA 

XA 

vO 

t^ 

o 

CM 

-d- 

C7^ 

r- 

^ 

vO 

C^ 

r- 

•LA 

lA 

o 

r^ 

O 

co 

lA 

(T- 

vO 

vO 

O 

CO 

CJ^ 

'- 

o 

CA 

\A 

t— 

\ 

° 

CM 

^ 

-•* 

• 

• 

• 

O 

CM 

lA 

-d- 

CO 

CO 

^ 

o 

CT^ 

-d- 

CO 

oo 

^ 

XA 

Cvj 

lA 

O 

-rt 

o 

vO 

vX) 

o 

^ 

o 

XA 

CM 

O 

•LA 

fA 

C7^ 

-d- 

CM 

CM 

-d- 

o> 

fA 

XA 

O^ 

C^ 

CO 

O 

O 

^— 

O 

O 

T— 

o 

o^^ 

CO 

t^ 

• 

« 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

- — ^ 

(M 

^— 

o 

r- 

o 

o 

o 

o 

^— 

o 

,_ 

CM 

vO 

'•-> 

■•-5 

•""5 

•<-> 

'f-J 

•"-i 

•f-s 

••~3 

••-3 

•<~i 

••-3 

•■"3 

t— 

+ 

+ 

1 

+ 

+ 

+ 

1 

1 

+ 

1 

1 

• — ' 

CM 

T— 

fA 

o 

CT^ 

-d- 

-d- 

c^ 

o 

fA 

CM 

O 

o 

lA 

C7^ 

CO 

V 

CO 

o 

VA 

o 

O 

O 

lA 

CA 

O 

r- 

5 

vD 

r- 

C7^ 

fA 

XA 

o 

o 

u\ 

-d- 

C^ 

CM 

r- 

\A 

C7^      ^      cr^ 

tA 

r--      CM      f- 

^ 

XA 

o 

CD 

o 

O 

CM 

o 

CD 

• 

o 

•                •                • 

O        CM        O 

• 

o 

• 

o 

' 

1 

1 

• 

1 

1 

1        1 

1 

vO 

lA 

CO 

CO 

m 

CO 

r- 

u^ 

XA 

[^ 

CM 

r-- 

CM 

CM 

r- 

CM 

f- 

o 

t— 

o 

r- 

r^ 

lA 

r- 

,— 

\A 

fA 

o 

nO 

o 

o 

vD 

o^ 

O 

CO 

o 

o 

00 

O 

C> 

vO 

o 

o 

O 

o 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

♦ 

^~^ 

O 

O 

r- 

r— 

o 

O 

r- 

o 

o 

o 

o 

,_ 

s 

••-i 

•r-S 

••-3 

'*": 

•f-s 

•'-i 

••-5 

•r-i 

•>-» 

••-3 

•"-3 

•r-i 

CO 

+ 

+ 

1 

+ 

1 

1 

+ 

+ 

■f 

1 

1 

1 

— 

CO 

\A 

sO 

o 

m 

CO 

r- 

CM 

r— 

CM 

r- 

o 

J- 

r- 

r- 

-^ 

o 

vO 

5 

CO 

r~ 

sD 

r^ 

T— 

vD 

v£> 

^— 

fA 

fA 

lA 

CO 

CO 

XA 

fA 

CM 

CM 

lA 

r- 

CM 

t^ 

LA 

CM        CM         O 

r~ 

o      o      o 

r-- 

O 

fvj 

^ 

O 

<D 

CM 

o 

O 

T-        CM        O 

O 

•                 •                  • 

o      o      o 

• 

o 

• 

o 

1 

• 

1 

' 

t 

vD 

CO 

CO 

CO 

r-- 

XA 

•LA 

r~- 

CM 

CM 

CM 

CM 

r- 

*— 

r- 

r^ 

r^ 

^— 

T- 

c<~\ 

o 

o 

o 

cr- 

CT^ 

o 

o 

o 

o 

O 

c^ 

c^ 

o 

o 

o 

a^ 

O 

^ 

^— 

r^ 

T~ 

O 

^ 

CM 

o 

• 

o 

CM 

• 

^~v 

•r-J 

•1-3 

•r-) 

■r-} 

•"-J 

••—> 

•r-s 

•r-} 

••-5 

•"-J 

•r-} 

•'-3 

CO 

+ 

+ 

1 

+ 

1 

1 

+ 

+ 

+ 

1 

1 

1 

' — ■ 

CO 

vD 

vO 

CO 

r>- 

<— 

r— 

r- 

o 

t^ 

r- 

o 

sO 

CO 

CO 

vD 

r^ 

v£) 

sO 

fA 

fA 

CO 

CO 

fA 

CM 

CM 

vO 

r- 

CM 

r- 

vD 

CM        CM         O 

fA 

o      o      o 

fA 

o 

CM 

r— 

o 

o 

1 

CM 
1 

o 
1 

O 

r-         CM         O 

1 

rA 

•                  •                 • 

o      o      o 

• 

fA 

• 

o 

CO 

CO 

CO 

CO 

^~ 

<r" 

o 

o 

CM 

CM 

J- 

-d- 

CO 

CO 

m 

■LA 

t^ 

r- 

(A 

fA 

o 

cr- 

CM 

CM 

CJ^ 

cr- 

C^ 

c^ 

C7^ 

o^ 

t^ 

r- 

• 

• 

• 

• 

• 

• 

* 

• 

'—» 

o 

o 

o 

o 

o 

o 

o 

o 

s 

••-J 

•■-3 

•'-J 

•'-5 

•1—5 

•r-> 

•l-J 

•'-3 

+ 

1 

1 

+ 

+ 

1 

+ 

1 

— 

sD 

vO 

rA 

fA 

-d 

-d 

r^ 

C^ 

<— 

T— 

CO 

CO 

CM 

CM 

CO 

CO 

ro 

<A 

r- 

r- 

v£) 

vO 

CO 

CO 

O 

CM 

vO 

CM 

CM 

vO 

O 

^X)         r-         O 

fA 

O        T-        vO 

Ch 

vJD 

o 

O 

o 
1 

CD 

CM 

o 

-- 

O        V-        o 

1 

fA 

o 

•                    • 

o 

CM 

O 

rH 
0) 
> 

0 

0) 

o 

o 

,— 

CM 

(A 

-d- 

■LA 

v£> 

t--        CC 

cr^ 

O 

T-               <\ 

r^. 

-d 

XA 

►-> 

u 

t— 

r-          r- 

r— 

T— 

^  ., 

4-> 
•H 

O 

bOvD 

r-^     I 

<  r- 


:3 

•H 

Cm 
O 

s 


32 


FFT         FFT  Normalized  to  mean  of  1.0 

D.C  0.0  0.0 

Tq  6.i^56i+  1.i;371 

2fQ  5.3933  1.2005 

3fQ  i|.3603  0.9705 

k^Q  2.2  O.I1897 

5fQ  3.i^.935  0.7776 

6fQ  3.1+973  0.7785 
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1.3559 
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^^0 
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0.14-756 

-^^0 

2.169/^- 

1  . 1  i;91+ 

Comparison  of  the  output  coefficients  of  the  FFT 
and  Modified  FFT  before  and  after  they  have  been 
given  a  mean  value  of  1 

Figure  1-7 
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shows  the  modified  FFT  algorithm  starting  from  the  4  level 
(the  FFT  and  modified  FFT  algorithms  are  the  same  up  through 
the  4  level).   Figure  1-7  compares  the  output  coefficients 
of  the  modified  FFT  and  the  FFT  before  and  after  the  coeffi- 
cients have  been  given  a  mean  value  of  1.   The  FFT  coeffi- 
cients have  also  been  given  a  mean  value  of  1  for  comparison 
purposes  which  shows  the  value  of  this  procedure  for  com- 
paring sequences  of  numbers  when  there  is  no  known  relation- 
ship between  the  two  sequences.   In  the  case  of  the  FFT  the 
mean  value  of  the  coefficients  was  computed  using  all  of  the 
coefficients  as  opposed  to  the  manner  in  which  it  was  per- 
formed in  the  case  of  the  modified  FFT.   There  is  no  loss  in 
generality  in  doing  this  because  dividing  a  sequence  by  its 
mean  value  in  no  way  alters  the  relative  magnitudes  of  the 
coefficients. 

It  is  interesting  to  note  in  the  case  of  the  modified 
FFT  that  the  signal  frequence  fo  (the  1  hertz  sinusoid)  ap- 
pears much  larger  than  2fo  before  the  "mean  1"  process  was 
applied.   After  applying  this  process  the  magnitude  of  fo  and 
2fo  are  very  close,  as  expected.   Grouping  the  output  coeffi- 
cients of  the  modified  FFT  that  depend  on  normalization  at 
each  level  and  then  comparing  them  with  the  FFT  coefficients 
it  is  seen  that  the  signals  of  the  modified  FFT  have  greater 
magnitude  than  the  signals  of  the  FFT  and  the  noise  coeffi- 
cients are  smaller  for  the  modified  FFT  than  for  the  FFT  in 
every  case  except  4fo. 
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Mod  FFT 

FFT 

Mod  FFT 

1.5031 

2fo 

1.2005 

< 

1.5243 

.6291 

6fo 

.7785 

> 

.4756 

FFT 

fo  1.4371  < 

3fo  .9705  > 

5fo  .7776  >    .7184 

7fo  1.3455  >   1.1494      4fo     .4897   <   1.0 

The  FFT  spectrum  was  given  a  mean  of  1  by  dividing  each 
coefficient  by  the  mean  of  all  of  the  coefficients.   In  the 
case  of  the  modified  FFT  each  coefficient  was  divided  by  the 
mean  of  the  coefficients  that  were  normalized  at  its  normal- 
ization level.   Each  procedure  resulted  in  a  spectrum  with  a 
mean  of  1  but   by  different  methods.   To  alleviate  any  doubt 
that  the  different  methods  used  to  give  each  process  a  mean 
value  of  1  were  responsible  for  the  apparent  improved  process 
gain  of  the  modified  FFT  over  the  FFT  in  this  example  the 
following  compares  the  processes  where  each  output  spectrum 
has  been  given  a  mean  value  of  1  by  the  method  used  for  the 
modified  FFT. 

FFT  Mod  FFT             FFT       Mod  FFT 

fo     1.2667  <  1.5031  2fo    1.2133  <   1.5243 

3fo     .8619  >  .6291  6fo     .7867  >    .4756 

5fo     .6854  <  .7184 

7fo    1.1860  >  1.1494  4fo    1.0  =   1.0 

The  method  of  normalizing  coefficients  to  a  mean  of  1 
seems  to  have  its  merits  but  still  does  not  solve  the  prob- 
lem of  giving  each  frequency  the  same  dependence  on  phase  co- 
herence as  opposed  to  amplitude  plus  phase  coherence  in  the 
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case  of  the  standard  FFT,   Consider  the  frequency  (N/4)fo. 
It  is  in  fact  the  FFT  coefficient  with  a  magnitude  of  1. 
The  frequencies  (N/S)fo  and  (3N/8)fo  do  not  have  their  ampli- 
tude influences  removed  until  just  prior  to  the  last  pass 
through  the  algorithm.   On  the  other  hand  the  frequencies 
nfo  (n=l , 3 , 5 , . . . , N/2-1)  have  their  amplitude  influences  re- 
moved from  their  intermediate  coefficients  early  in  the 
transform. 

D.   METHOD  OF  IMPLEMENTATION 

The  method  that  was  finally  chosen  was  that  of  making  all 
frequencies  look  like  odd  harmonics  of  1/T.'  Keeping  T  con- 
stant and  using  frequency  translation  was  tried  but  did  not 
solve  the  problem  as  the  even  harmonics  of  fo,  when  trans- 
lated by  1/T,  still  behaved  like  even  harmonics  of  fo  when 
normalization  was  performed  at  the  4  level.   That  is,  they 
did  not  reflect  all  possible  phase  angles  at  the  output  of 
the  transformation. 

The  method  that  did  work  was  to  vary  T,  the  length  of  the 
transform,  so  that  all  frequencies  were  made  to  look  like  odd 
harmonics  of  1/T.   This  is  best  explained  by  example.   Con- 
sider a  time  record  T  seconds  long  and  a  sinusoid  of  fre- 
quency 2/T  that  has  been  sampled  N  times  in  T  seconds. 

The  correlating  sinusoids  are: 
cos(2TTnk/N)-jsin(2TTnk/N) . 
The  only  non-zero  correlation  occurs  when  n=2.   Now,  if  the 
time  record  were  divided  in  half,  then  there  would  be  N/2 
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samples  in  time  T/2.   The  correlating  sinusoids  would  then 
be:  cos[2TTnk/(N/2)]-jsin[2Trnk/(N/2)]  and  the  only  non-zero 
correlation  would  occur  when  n=l.   The  net  result  is  that 
the  correlating  sinusoids  are  identical  but  the  transform 
length  is  cut  in  half  and  the  2fo  hertz  sinusoid  looks  like 
a  sinusoid  of  frequency  fo  to  the  time  record  of  length  T/2. 
If  two  transformations  were  performed  and  the  output  coeffi- 
cients are  summed  in  the  complex  domain  the  net  result  is 
the  same  as  if  one  transformation  was  performed  of  length  T. 
This  is  simply  because  the  correlating  sinusoids  were  identi- 
cal in  both  cases,  that  is,  every  multiplication  was  identi- 
cal in  each  case  as  was  the  summation  of  the  products  of 
these  multiplications.   If  a  4fo  sinusoid  were  made  to  look 
like  an  odd  harmonic,  then  the  input  record  would  have  to  be 
divided  into  four  equal  parts  and  4-N/4  point  transformations 
would  then  have  to  be  summed. 

When  the  normalization  modification  is  then  added  to  the 
algorithm  essentially  only  one-half  of  the  algorithm  is  per- 
formed, the  half  that  deals  with  the  odd  harmonics  of  fo. 
But  the  algorithm  must  be  performed  considerably  more  times. 
First  all  N  points  are  transformed  to  obtain  the  frequency 
coefficients  relating  to  nfo  (n=l , 3 , 5 , .  .  .  , N-1 ) .   Then  2-N/2 
point  transforms  are  performed  to  obtain  the  frequency  coef- 
ficients relating  to  nfo  (n=2 , 6 , 10, . . . , N-2) .   Then  4-N/4 
point  transforms  are  performed  to  obtain  the  frequency  coef- 
ficients relating  to  nfo  (n=4 , 12 , 20, . . . , N-4 ) ,  etc.,  until 
N/4-4  point  transforms  are  performed  to  obtain  the  frequency 
coefficient  relati;u;  to  (N/4)fo. 
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For  a  1024  point  transform  the  time  required  to  obtain 
all  of  the  frequency  coefficients  is  approximately  three 
times  that  of  a  normal  FFT  if  the  transformations  are  per- 
formed one  after  the  other.   If  parallel  processing  were  em- 
ployed the  processing  time  could  be  reduced  to  half  of  that 
required  for  the  FFT.   If  the  criterion  of  NLogzN  (for  the 
nvmiber  of  complex  operations  required  in  a  transformation)  is 
used  for  comparing  the  relative  time  required  for  the  two 
transformations  then,  for  32  points  it  requires  32  log2(32)  = 
32x5  =  160  for  the  FFT.   For  the  modified  FFT  it  takes 
16x(5+4+3+2)  =  16x14  =  224.   For  a  1024  point  transform  the 
FFT  requires  1024x10  =  10,240  where  the  modified  FFT  re- 
quires 512x(10+9+. . .+2)  =  512x54  =  27,648,   For  8192  points 
the  FFT  requires  106,496  and  the  modified  FFT  requires 
368,640.   The  numbers  do  not  account  for  the  extra  bit  re- 
versals, additional  ordering  of  output  coefficients,  or  the 
normalization  procedure  required  in  the  modified  FFT. 

To  summarize  the  results  of  this  transformation:   The 
proposed  modification  to  the  FFT  algorithm  will  be  more  de- 
pendent on  the  phase  coherence  of  a  waveform  to  be  trans- 
formed and  less  dependent  on  the  power  in  the  waveform  than 
the  FFT.   When  a  sinusoid  of  the  form  Asin(2TTnf  ot  +  0 )  (n=l,2, 
3,..., N/2-1)  is  transformed  the  resulting  output  coefficients 
will  be  a  constant  independent  of  A,  n,  or  0.   And  9  can  be 
recovered  from  the  output  coefficient  corresponding  to  n. 

The  performance  of  this  algorithm  has  been  thoroughly  in- 
vestigated for  signals  of  the  form  in  the  preceding  paragraph, 
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Due  to  the  nonlinearity  of  the  process  it  is  difficult  to 
predict  how  it  will  perform  in  general.   The  next  section  in- 
vestigates the  performance  of  the  algorithm  against  that  of 
the  FFT  under  identical  conditions  for  different  types  of 
data. 
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III.   PART  2  EXPERIMENTAL  RESULTS 

A.   MEASUREMENTS 

In  determining  the  value  of  the  modification  to  the  FFT 
various  types  of  data  were  transformed  and  compared  v/ith  the 
FFT  transformations.   Several  measurements  were  employed  to 
evaluate  the  relative  merits  of  the  FFT  modification  when 
compared  to  the  FFT.   The  output  spectrum  consisted  of  the 
magnitudes  of  the  output  coefficients  referred  to  as  {F.}, 
(i=0,l,2, . . ,N/2-l),   A  coefficient  that  is  related  to  the 
frequency  of  an  input  signal  is  referred  to  as  F  .   In  gen- 
eral,  the  objective  of  signal  processing  is  the  identifica- 
tion of  a  coefficient  that  is  signal  related  from  one  that 
is  not  signal  related.   Norm.ally  this  is  accomplished  by 
visual  observation  of  the  magnitude  of  the  output  spectrums 
of  the  transformations.   Large  coefficients  tend  to  indicate 
the  presence  of  a  signal  and  just  as  important,  the  consist- 
ency of  the  amplitude  of  an  output  coefficient  from  trans- 
form to  transform  indicates  the  presence  of  a  signal.   It  is 
the  combination  of  these  two  effects  that  gives  reliability 
to  the  estimate  of  signal  presence  or  absence. 

What  may  be  one  of  the  best  means  of  determining  the 
presence  or  absence  of  a  signal  is  the  analysis  by  an  exper- 
ienced spectrum  analyst.   Amplitude  and  consistency  of  an 
output  coefficient  do  not  always  portray  all  of  the  informa- 
tion available  in  a  coefficient.   There  seem  to  be  other 
characteristics  in  the  behavior  of  a  spectrum  from  transform 
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to  transform  that  gives  extra  support  to  an  analyst's  abil- 
ity to  determine  if  a  signal  is  present  or  absent.   In  gen- 
eral it  is  difficult  to  explain  what  these  additional 
characteristics  are  and  it  seems  that  only  after  a  great 
deal  of  experience  do  they  figure  in  to  an  analyst's  logic 
when  he  decides  if  a  signal  is  present  or  not.   The  extra 
signal  characteristics  generally  pertain  to  a  specific  type 
of  signal  in  a  particular  environment  with  which  the  analyst 
is  very  familiar;  that  is,  he  cannot  generalize  his  experi- 
ence with  a  specific  class  of  signals  to  all  types  of  sig- 
nals and  environments. 

The  primary  methods  utilized  in  determining  the  merits 
of  the  modified  FFT  were: 

1.  Estimation  of  signal  enhancement  through  visual  ob- 
servation . 

2.  Determining  the  difference  in  amplitude  between  the 
magnitude  of  the  signal  and  the  magnitude  of  the  mean  value 
of  the  noise: 

F i^{[  E  F.  ]-F  }  =  F   -  p 

s    n-1^  ^ .^^  1^   s'     s    ^n 

where  u   is  the  mean  value  of  the  noise  coefficients, 
^n 

3.  Determining  the  variance  of  the  noise: 

normally  associating  the  smaller  variance  with  the  better 
processing  gain  (referred  to  as  smoothing  the  noise). 
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4.  Utilization  of  a  threshold  detector  to  determine 

with  what  probability  F  would  exceed  all  of  the  other  F.  or 

all  of  the  F.  within  a  neighborhood  of  F  ,  such  as  F  r.<r .  • 
1  "^  s  s-25+1 

(i=0, 1,2, . . . ,50) ,  where  the  neighborhood  in  this  case  is  25 

coefficients  on  either  side  of  F  . 

s 

5.  Comparing: 

(F  -  a  )2 
s    n 


2a2 
n 

for  the  different  spectrums.   The  numerator  relates  how  much 
greater  F   is  than  the  average  of  the  comparison  noise  coef- 
ficients  and  the  denominator  relates  to  the  smoothness  of 
the  noise  coefficients.   The  form  of  this  equation  is  identi- 
cal to  that  normally  used  in  determining  input  signal  to 
noise  ratios. 

Method  5  was  used  in  almost  every  experiment  (it  was  not 
used  in  some  of  the  early  experiments)  along  with  one  or 
more  of  the  other  methods.   It  was  found  that  methods  4  and 
5  agreed  quite  closely  in  almost  every  case  (there  were  a 
negligible  number  of  exceptions  and  then  the  difference  it- 
self was  negligible).   "Agreed"  in  the  context  of  the  above 
sentence  means  that  if  a  transformation  was  better  than 
another  in  one  test  it  was  also  better  in  the  other  test. 
Method  5  is  the  easiest  to  compute  and  relates  directly  to 
the  signal  to  noise  ratio  of  the  output  spectrum.   As  a  re- 
sult it  is  used  most  often  in  the  following  pages  to  make 
comparisons  of  the  FFT  and  the  modified  FFT. 
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The  problem  of  determining  input  signal  to  noise  ratios 
across  a  particular  noise  bandwidth  was  alleviated  by  com- 
paring the  output  spectrum  of  the  modified  FFT  to  that  of 
the  FFT.   In  this  way  the  performance  of  the  modification 
was  judged  relative  to  the  FFT;  i.e.,  the  FFT  was  used  as 
the  standard  of  comparison. 

It  was  originally  hoped  that  the  resulting  spectrum  from 
the  modified  FFT  would  differ  significantly  from  that  of  the 
standard  FFT.   If  that  was  true,  and  if  the  modified  FFT  had 
sufficient  process  gain  then  the  combination  of  the  two 
spectra  could  result  in  a  spectrum  with  a  lower  a   and  then 
hopefully  a  greater  process  gain. 

Since  all  of  the  theoretical  work  presented  so  far  deals 
with  simple  sinusoids  at  multiples  of  the  fundamental  fre- 
quency of  the  transformation  a  large  void  is  left  between 
what  happens  in  that  case  and  in  the  case  of  an  input  wave- 
form consisting  of  multiple  sinusoids  at  arbitrary  frequen- 
cies combined  with  random  noise,  Gaussian  or  not,  especially 
since  the  process  is  non-linear. 

Mathematical  or  statistical  analysis  of  the  modified  FFT 
algorithm  is  complicated  at  best,  and  since  the  process  is 
non-linear  all  conclusions  as  to  its  value  were  developed 
through  experimentation  with  different  types  of  input  data. 
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B.   SINGLE  SINUSOID  PLUS  GAUSSIAN  NOISE 

The  first  data  chosen  for  analysis  consisted  of  a  simple 
sinusoid  combined  with  additive  Gaussian  noise.   Most  noise 
generators  are  designed  to  generate  Gaussian  noise  and  its 
statistical  properties  are  well  known  making  it  the  most  con- 
venient and  generally  accepted  type  of  data  for  comparison 
purposes.   Due  to  the  simplicity  of  a  single  sinusoidal  sig- 
nal and  the  availability  of  Gaussian  noise  generators  the 
first  experiments  were  performed  on  this  type  of  data. 

The  signal  chosen  was  a  25  hz  sinusoid  with  additive 
Gaussian  noise  with  a  100  hz  bandwidth.   The  signal  plus 
noise  was  generated  with  a  "Wavetek"  model  132  VCG/Noise  Gen- 
erator.  The  input  waveform  was  lowpass  filtered  at  100  hz 
and  sampled  at  256  samples  per  second.   The  transformations 
were  2048  points  in  length  (8  seconds  in  time)  with  a  re- 
sulting resolution  of  1/8  hz.   The  noise  to  signal  ratios 
ranged  from  0  to  24  db  in  3  db  increments,  for  3.33  minutes 
or  25,  8-second  transforms  at  each  level. 

Figures  2-1  through  2-9  show  the  results  of  these  trans- 
formations.  The  horizontal  axis  represents  frequency,  in- 
creasing from  left  to  right.   150  points  (output  coefficients) 
are  displayed.   The  center  of  the  highlighted  area  represents 
the  25  hz  frequency.   The  left  edge  is  15.625  hz  and  the 
right  edge  is  34.375  hz.   The  vertical  axis  represents  the 
relative  magnitude  of  the  coefficients.   The  fore  and  aft 
axis  represents  time  with  the  most  current  information  being 
displayed  at  the  front  (or  bottom)  and  the  oldest  information 
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being  displayed  at  the  back  (or  top).   Each  trace  in  each 
figure  represents  the  average  of  ten  transforms  or  80  seconds 
of  data.   Each  trace  contains  80%  of  the  same  information 
contained  in  the  traces  adjacent  to  it.   A  boxcar  type  aver- 
ager is  used  with  five  compartments.   In  this  case  two  trans- 
forms were  averaged  in  each  compartment,  then  the  five 
compartments  were  averaged.   Each  succeeding  trace  is  formed 
by  dumping  the  compartment  with  the  oldest  information  (20% 
of  the  information  in  the  boxcar)  and  refilling  it  with  new 
information,  averaging  and  then  displaying.   With  this  in 
mind  then  only  traces  that  are  separated  by  five  other  traces 
are  independent  in  terms  of  time  from  one  another. 

Table  I  gives  the  signal  to  noise  output  ratios  computed 
for  the  average  of  25  transforms  at  each  db  level.   The  noise 
statistics  were  computed  from  the  25  coefficients  on  each 
side  of  the  signal  coefficient.   The  coefficients  immediately 
to  the  left  and  right  of  the  signal  were  not  included  in 
these  computations.   As  stated  earlier  all  S/N  output  ratios 
were  computed  by  subtracting  the  mean  of  the  noise  coeffi- 
cients from  the  signal  coefficient  and  dividing  the  square 
of  that  result  by  twice  the  variance  of  the  noise  coeffi- 
cients. 

The  most  obvious  result  of  comparing  the  respective  out- 
put spectra  of  the  FFT  and  modified  FFT  is  their  strong  simi- 
larities.  Close  observation  of  the  figures  reveals  the 
almost  identical  behavior  of  the  respective  spectra  as  the 
S/N  ratio  is  reduced.   Also  the  relative  magnitude  of  the 


45 


signal  coefficient  with  respect  to  the  magnitude  of  the  co- 
efficients immediately  adjacent  to  the  signal  coefficient 
seems  to  be  the  same  for  both  spectra.   This  effect  adds 
another  element  of  similarity,  the  effect  of  having  the  same 
"shape".   The  observed  variance  of  the  noise  appears  to  be 
somewhat  greater  in  the  modified  FFT  spectrum. 

Table  I  readily  confirms  the  observed  similarities  in 
the  spectra.   The  S/N  outputs  for  the  two  spectra  are  very 
close.   The  FFT  generally  has  better  S/N  ratios  and  this  can 
be  attributed  to  the  smaller  noise  variance  observed  in  the 
photographs.   The  FFT  becomes  significantly  better  than  the 
modified  FFT  for  the  -18,  -21,  and  -24  db  S/N  ratios,  this 
is  again  confirmed  in  the  photographs.   Comparison  of  the 
FFT  spectrum  squared  with  that  of  the  FFT  spectrum  multiplied 
by  the  modified  FFT  spectrum  shows  no  signal  enhancement 
over  what  can  be  gained  by  simply  squaring  the  FFT  spectrum. 
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FFT 


Modified  FFT 

25  hz  signal,  S/N  =  0  db 

10  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-1 
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Modified  FFT 

25  hz  signal,  S/N  =  -3  db 

10  transformed  averaged  per  trace,  transform  length  204  8 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-2 
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Modified  FFT 


25  hz  signal,  S/N  =  -6  db 

10  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 


Figure  2-3 
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FFT 


Modified  FFT 

25  hz  signal,  S/N  =  -9  db 

10  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-4 
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FFT 


Modified  FFT 

25  hz  signal,  S/N  =  -12  db 

10  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-5 
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FFT 


Modified  FFT 

25  hz  signal,  S/N  =  -15  db 

10  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-6 
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FFT 


Modified  FFT 

25  hz   signal,  S/N  =  -18  db 

10  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-7 
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FFT 


Modified  FFT 

25  hz  signal,  S/N  =  -21  db 

10  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-8 


54 


FFT 


Modified  FFT 

25  hz  signal,  S/N  =  -24  db 

10  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-9 
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S/N  out  (db) 

^/!Jk^"       FFT      Mod-FFT  FFTxFFT  ./^"^Jx^m 

(db  Mod-FFT 

0      26.295     27.260  39.790  40.902 

-3      29.734     30.106  42.409  42.414 

-6      35.007     32.860  48.469  47.258 

-9       34.320     34.173  47.235  47.578 

-12      25.480     24.053  32.632  32.208 

-15      21.383     20.342  26.925  26.656 

-18      21.311     18.646  26.124  25.006 

-21       11.519      7.664  13.726  12.490 

-24       9.877      2.422  11.457  8.506 

Table  I 

Average  S/N  out  (db)  vs.  S/N  in  (db)  for  the  Average  of  25 

Transforms  in  the  case  of  the  25  hz  signal  plus  Gaussian 

noise.   S/N  out  computed  as: 

(F  -VI  )^ 
^  s   n 


2o^ 
n 

where  F   is  the  magnitude  of  the  signal  coefficient  and  y 
s  ^  "^^  n 

and  a   are  computed  from  25  noise  coefficients  on  either  side 

of  F  . 
s 
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C.   MULTIPLE  SINUSOIDS  PLUS  GAUSSIAN  NOISE 

To  further  investigate  the  relative  merits  of  the  modi- 
fied FFT  Gaussian  noise  was  again  used.   In  this  case  multi- 
ple sinusoidal  signals  were  present  allowing  investigation 
into  the  possibility  of  destructive  interference  between  the 
signals  due  to  the  normalization  procedure. 

The  data  came  from  a  test  tape  made  by  the  Naval  Air 
Development  Center.   The  resulting  spectra  are  presented  in 
Figures  2-10  through  2-14.   The  data  consisted  of  four  sig- 
nals in  two  groups  of  two.   Only  one  group  of  two  was  consid- 
ered.  The  tape  had  ten  minutes  of  data  for  each  signal 
starting  with  three  minutes  of  white  noise  followed  by  ten 
minutes  of  signals  at  -20  db  and  -18  db  presented  in  close 
proximity  to  each  other.   These  were  followed  by  two  more 
signals  for  ten  minutes  at  -16  db  and  -14  db  etc.,  until  two 
signals  were  presented  at  0  db  and  +2  db  (the  last  two  were 
not  transformed).   The  input  signal  to  noise  ratios  were 
measured  with  respect  to  a  1  hz  noise  bandwidth.   The  fig- 
ures indicate  where  the  signals  are  or  where  they  should  be. 

These  figures  are  interpreted  in  exactly  the  same  manner 
as  Figures  2-1  through  2-9.   In  this  case  forty  transforms 
of  2048  points  each  were  averaged  in  each  trace.   The  boxcar 
average  was  again  used.   This  data  was  also  filtered  at  100 
hz  and  sampled  at  256  samples  per  second.   The  frequency 
resolution  is  1/8  hz. 

It  should  be  noted  that  there  was  some  question  as  to 
the  reliability  of  the  -18  db  signal  and  it  was  not  considered 
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to  be  experimentally  acceptable  in  the  data  analysis.  The 
indication  of  where  it  should  be  in  Figures  2-14  is  only  a 
best  guess. 

Table  II  gives  the  average  computed  output  S/N  ratios 
for  the  summation  of  40  transforms.   These  were  computed 
using  the  same  method  employed  throughout  this  study.   In 
this  case,  when  a  signal  was  on  the  right  only  25  coeffi- 
cients to  the  right  of  that  signal  were  considered  in  com- 
puting the  noise  statistics.   Likewise,  if  the  signal  was  on 
the  left  only  25  coefficients  to  the  left  of  that  signal 
were  considered  in  computing  the  noise  statistics.   In  all 
cases  the  nearest  coefficient  to  the  signal  coefficient  was 
not  considered  in  the  computations. 

Table  III  is  similar  to  Table  II. except  that  80  trans- 
forms were  averaged  instead  of  40. 

The  most  important  result  of  this  experiment  was  to  show 
that  multiple  signals  did  not  interfere  with  each  other  and 
that  the  normalization  process  did  not  appear  to  destroy  or 
reduce  signal  resolutions.   Comparison  of  the  figures  again 
reveals  similar  spectra  for  the  respective  transformations. 
Comments  identical  to  those  made  in  Section  III-B  with  re- 
spect to  the  single  sinusoid  plus  Gaussian  noise  case  could 
be  made  here  with  one  exception.   The  spectra  for  the  -16 
and  -20  db  signal  to  noise  ratios  tend  to  favor  the  modified 
FFT.   In  Figures  2-13  the  top  trace  of  the  modified  FFT 
shows  the  -16  db  signal  to  be  nearly  as  large  as  the  -14  db 
signal  and  in  the  second  trace  from  the  top  it  is  actually 
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larger.   Also,  the  variance  of  the  adjacent  noise  coeffi- 
cients appears  to  be  smaller  in  the  modified  FFT  spectrum 
than  in  the  FFT  spectrum.   The  individual  distinctions  are 
not  as  pronounced  in  the  case  of  the  -20  db  signal  in  Fig- 
ures 2-14  but  the  general  appearance  tends  to  favor  the  mod- 
ified FFT  spectrum. 

Table  II  indicates  that  the  FFT  generally  outperformed 
the  modified  FFT  when  the  signal  to  noise  ratios  were  rela- 
tively small.   The  reverse  was  true  for  the  -16  and  -20  db 
S/N  ratios  (as  indicated  above  the  -18  db  signal  is  not  being 
considered  in  the  analysis).   As  in  the  case  of  the  single 
sinusoid  presented  in  Section  III-B  no  obvious  trends  were 
present  in  this  data.   Again  the  product  of  the  two  spectra 
when  compared  against  the  FFT  spectrum  squared  did  not  prove 
to  have  any  particular  advantage. 

Table  III  compares  the  measured  results  when  80  trans- 
forms were  averaged  instead  of  40  as  in  Table  II.   The  im- 
provement in  S/N  output  was  relatively  consistent  for  all  S/N 
ratios  and  for  both  spectra  with  neither  spectrum  gaining  any 
particular  advantage  over  the  other  as  a  result  of  the  deeper 
average. 
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I|.0   transforms    averaged  per   trace,    transform  length  20i|.8 
Displayed   bandwidth   18,75  hz,    0.12^  hz   resolution 
.   -2  db   signal   on   right,    -l^.  db   signal   on   left 


H^DC   TlST   TflPE   (rPTj 


S/N--2, -4  DB 


Figure   2-10    (FFT) 
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L|.0  transforms  averaged  per  trace,  transform  length  20J|8 
Displayed  bandwidth  18.75  hz,  0,123'  hz  resolution 
-2  db  signal  on  right,  -l\.   db  signal  on  left 


X-jf.ALc   -    -Xh^-jl-^i  LiriJTj/JNCM 


^-oC.mLc   -  3..m^'fil  U^JTo'JNCH, 


imc  TfST  ta^e;  (mod  fptj 

3/H^  ~2     "4     DB  Figure    2-10    (Kod-FFT) 
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1|0  transforms  averaged  per  trace,  transform  length  20[t.8 
Displayed  bandwidth  18.75  hz,  0,125  hz  resolution 
-6  db  signal  on  left,  -8  db  signal  on  right 


NqDC   TTST   Tfl^-^f:   (rPTj 
S/N---6,  -8  DB  f'ig^^^  2-11  (FFT) 
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1^.0   transforms  averaged  per  trace,  transform  length  20b,8 
Displayed  bandwidth  18.75  hz,  0.12$  hz  resolution 
-6  db  signal  on  left,  -8  db  signal  on  right 


«C   TTST   TfiPL  (MOD  FPT) 
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w^wnl 


[|.0  transforms  averaged  per  trace,  transform  length  20I4.8 
Displayed  bandwidth  18.75  hz,  0,125  hz  resolution 
-10  db  signal  on  left,  -12  db  signal  on  right 


mx  TfST   TflPE  (PFTj 

S/N---10P  -12     DB         Figure    2-12    (FPT) 
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I4.O  transforms  averaged  per  trace,  transform  length  20)4.8 
Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 
-10  db  signal  on  left,  -12  db  signal  on  right 


X-:)C.^LI  -  Jmi-^l   LJ.^JTj/jriC!i, 

N^GC   TF;ST   Tfl^r   (MOD  FfTj- 

S/h--  10,-12     OB  Figure    2-12    (Mod-FFT) 
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I4.O  transforms  averaged  per  trace,  transform  length  20i4.8 
Displayed  bandwidth  18,75  hz,  0.12^  hz  resolution 
-II4;  db  signal  on  left,  -16  db  signal  on  right 


iqDC   TF^ST   TflDE   (PPTJ 

5/|\|-.  -  ]Lj^  -  jg     Qg  Figure   2-13    (FFT) 
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i|0  transforms  averaged  per  trace,  transform  length  20I4.8 
Displayed  bandwidth  18,75  hz,  0,12^?  hz  resolution 
-II4.  db  signal  on  left,  -16  db  signal  on  right 


^-5MLE   -  3.MZ-Pi]  U.'NJTj/JNCM, 

M^GC   TEST   Tflt^F.   (MOD  FPTj- 

5/[\|-..  -  ]LJ^  -  Jg     QQ  Figure   2-13    (Mod-FFT) 
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I4.O  transforms  averaged  per  trace,  transform  length  20i|8 
Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 
-18  db  signal  on  right,  -20  db  signal  on  left 


i^DC   TTST   TflDf   (pfTj 

3/N--I8,   ~20     DB  Pig^^e    2-11+    (FFT) 
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l;0  transforms  averaged  per  trace,  transform  length  20I4.8 
Displayed  bandwidth  18.75  hz,  0,125  hz  resolution 
-18  db  signal  on  right,  -20  db  signal  on  left 


X-:.L"^LE    -   3.MZ-^1  UNJT5/JHCH, 

^-5rflLE  -  3mz-!di  ij.n:to'Jmch. 

mC  TfST   Tfi^E  (MOD  PrTj 

/N^'IS,   -20     DB         Figure    2-11+    (Mod-FFT) 
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S/N  out  (db) 

^/!Jv.n"       FFT      Mod-FFT  FFTxFFT  ./^"^^vrn 

(db)  Mod-FFT 

-2      19.655     18.059  23.985  23.165 

-4       17.086     16.388  20.218  20.092 

-6                  2.064             5.086  0.358  2.521 

-8                  5.901             8.397  4.622  7.018 

-10                  6.882             5.693  4.769  5.373 

-12                  8.215             7.121  6.761  6.641 

-14                  8.847             7.258  7.499  7.160 

-16                   5.489              6.386  4.008  '    5.015 

-18                  5.399             8.119  3.086  5.394 

-20                   4.370              6.872  3.546  5.865 

Table    II 

Average  S/N  out  (db)  vs.  S/N  in  (db)  for  the  Average  of  40 

Transforms  from  the  NADC  Test  Tape.  S/N  computed  as: 

(F  -p  )^ 
^  s   n 


2a^ 
n 

where  F   is  the  magnitude  of  the  signal  coefficient  and  y 
s  '^  ^  n 

and  o   are  computed  from  25  noise  coefficients  on  one  side 

of  F  . 
s 
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S/N  out  (db) 

S/N  m       pp„      Mod-FFT  FFTxFFT  FFTx 

(db)       ^^^     ^°^     ^  mxrta  Mod-FFT 

-2       23.669     22.083  28.316  27.606 

-4  .    19.943     19.738  23.094  23.300 

-6       3.334      7.066      1.217  3.991   • 

-8      10.492     10.054      8.871  9.427 

-10                   8.575              8.366              6.420  7.140 

-12      10.395      8.889      8.827  8.437 

-14      10.895     10.048     9.414  9.021 

-16       7.355      9.487      5.707  7.171 

-18        7.504     11.319      4.676  7.294 

-20       7.032     10.372      6.065  8.581 

Table  III 

Average  S/N  out  (db)  vs.  S/N  in  (db)  for  the  Average  of  80 

Transforms  from  the  NADC  Test  Tape.  S/N  computed  as: 

(F  -y  )2 
^  s   n 


2o' 
n 

where  F   is  the  magnitude  of  the  signal  coefficient  and  y 

and  o   are  computed  from  25  noise  coefficients  on  one  side 

of  F  . 
s 
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D.   SINGLE  SINUSOID  PLUS  GAUSSIAN  NOISE  PLUS  BURST  NOISE 

The  data  to  be  investigated  in  this  section  is  similar 
to  that  discussed  in  Section  III-B,  except  in  this  case  addi- 
tional noise  was  added  in  the  form  of  random  bursts  spaced 
at  random  intervals  with  random  duration.   The  data  was  gen- 
erated in  this  manner  to  accentuate  the  effect  of  spectrum 
dominance  by  vectors  with  large  amplitudes.   This  type  of 
data,  where  large  bursts  of  noise  could  dominate  or  corrupt 
the  resulting  spectrum,  provides  a  good  test  of  the  FFT  mod- 
ification.  If  the  modification  is  performing  as  anticipated 
it  will  reduce  the  spectrum  distortions  caused  by  the  large 
amplitudes  of  the  burst  noise  by  weighting  all  of  the  inter- 
mediate coefficients  equally.   Then  large  output  coefficients 
will  only  be  obtained  through  phase  coherent  vector  additions 
and  not  through  vectors  with  large  amplitudes. 

In  this  case  the  signal  was  set  at  10  hz  vice  25  hz  and 
additional  noise  was  added  in  the  form  of  bursts  with  random 
length  and  spaced  at  random  intervals.   Figures  2-15  through 
2-17  give  samples  of  the  input  waveform.   In  Figure  2-15  the 
top  trace  shows  the  signal  plus  the  burst  noise,  the  bottom 
figure  shows  it  after  it  has  been  filtered  at  100  hz.   Fig- 
ures 2-16  and  2-17  show  the  burst  noise  alone  in  the  top 
trace,  the  signal  plus  noise  plus  burst  noise  in  the  middle 
trace  and  the  bottom  trace  is  the  middle  trace  after  it  has 
been  filtered  at  100  hz.   In  each  figure  the  time  scale  is 
0.25  seconds  per  division. 
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As  was  previously  donfe  the  data  was  sampled  at  256  sam- 
ples per  second  and  2048  point  transforms  were  then  computed. 
Figures  2-18  through  2-21  show  the  results  of  the  transforma- 
tions.  Each  trace  contains  the  average  of  40  transforms  and 
as  was  done  previously  there  is  80%  redundancy  between  adja- 
cent traces.   The  S/N  ratios  on  the  figures  refer  to  the  S/N 
ratios  of  the  signal  to  continuous  noise  and  do  not  account 
for  the  burst  noise. 

Figures  2-22  through  2-28  are  photographs  of  the  same 
material.   The  photographs  are  presented  because  the  graph 
routine  that  drew  Figures  2-18  through  2-21  would  not  func- 
tion if  fewer  than  five  traces  were  to  be  drawn.   Figures 
2-22  through  2-24  present  the  transforms  in  the  same  manner 
as  Figures  2-18  through  2-21.   Each  trace  is  the  summation 
or  average  of  40  transforms.   Figures  2-25  through  2-28 
present  the  results  of  averaging  80  transforms.   Figure  2-25 
shows  the  -24  db  S/N  transforms  and  Figure  2-26  is  the  same 
presentation  only  in  this  case  the  coefficients  have  been 
squared.   It  is  felt  that  squaring  a  spectrum  gives  better 
insight  into  the  signal  recognition  problem.   If  the  signal 
is  readily  discernable  before  the  coefficients  are  squared, 
then  after  squaring  it  will  be  even  more  discernable.   If 
the  signal  becomes  less  discernable  after  squaring,  then  it 
is  concluded  that  the  S/N  ratio  is  insufficient  to  allow  any 
degree  of  signal  recognition. 

Figures  2-27  and  2-28  present  the  -30  db  transforms 
exactly  as  described  for  the  -24  db  transforms. 
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Tables  IV  and  V  give  the  signal  to  noise  output  ratios 
for  the  different  S/N  ratios.   Table  IV  is  the  average  S/N 
output  ratio  for  the  summation  of  40  transforms,  and  Table  V 
is  for  the  summation  of  80  transforms.   In  each  of  these 
cases  the  noise  statistics  were  computed  using  25  coeffi- 
cients on  either  side  of  the  signal  coefficient,  not  includ- 
ing the  coefficients  adjacent  to  the  signal  coefficient. 

The  analysis  of  the  transformations  presented  in  this 
section  lead  quite  conclusively  to  the  result  that  the  modi- 
fied FFT  is  capable  of  providing  better  S/N  output  ratios 
than  the  FFT  on  this  type  of  data.   The  modified  FFT  provides 
substantially  improved  S/N  output  ratios  over  the  FFT  in 
Figures  2-18,  2-19,  and  2-20.   In  Figures  2-21  neither  trans- 
formation provides  sufficient  signal  recognition. 

Figures  2-22  through  2-24  are  photographs  of  the  spectra 
presented  in  Figures  2-18  through  2-20.   The  photographs 
emphasize  the  extra  noise  smoothing  provided  by  the  modified 
FFT.   Figures  2-25  through  2-28  result  from  the  average  of 
80  transforms  instead  of  40  as  was  done  in  the  previous  fig- 
ures.  In  each  of  these  figures  the  variance  of  the  noise 
coefficients  is  considerably  greater  in  the  FFT  spectrum 
than  in  the  modified  FFT  spectrum.   In  each  of  the  figures 
the  signal  is  more  distinct  and  more  easily  recognized  in 
the  modified  FFT  spectrum  than  in  the  FFT  spectrum. 

Table  IV  shows  the  consistently  better  S/N  output  ratios 
of  the  modified  FFT  over  the  FFT.   The  -24  db  input  S/N 
ratio  is  the  only  exception.   Another  interesting  observation 
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is  the  improvement  of  the  S/N  output  ratios  obtained  by  mul- 
tiplying the  FFT  and  modified  FFT  spectra  together.   The 
product  of  the  two  spectra  is  consistently  better  than  either 
the  individual  spectra  or  the  FFT  spectrum  squared. 

Table  V  further  emphasizes  the  modified  FFT's  merits  in 
transforming  this  data.   By  averaging  80  transforms  the  mod- 
ified FFT  was  able  to  produce  better  S/N  output  ratios  than 
the  FFT  in  every  case  tested. 
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UO  trsmsforms  averaged  per  trace,  transform  length  2014.8 
Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 


1+8UR3T     FP'T 
i--r2  DB 


Figure   2-18    (PPT) 
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I4.O  transforms  averaged  per  trace,  tranaform  length  20L1.8 
Displayed  bandwidth  18.75  hz,  0.1 2^  hz  resolution 


S-fM-fBURST     MOD  FFT 

S-/ri""~i2     DB  Figure    2-18    (Kod-PFT) 
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kO   transforms  averaged  per  trace,  transform  length  20l|8 
Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 


S-+M-fBURST     FFT 

S/M"^  "lo     DB  Figure    2-19    (FFT) 
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llO   transforms  averaged  per  trace,  transform  length  201^.8 
Displayed  bandwidth  18.75  hz,  0,125  hz  resolution 

K'5{.-^LZ  =  j,.F)f)i:-?i^  mjT.j/irscH. 


S/M--18  DB 


Figure    2-19    (Mod-PFT) 
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ij.0  transforms  averaged  per  trace,  transform  length  20I4.8 
Displayed  bandwidth  18,75  hz,  0,125  hz  resolution 


-+M-+BU^5T     FPT 


lO 


S--h--24  Db 


Figure   2-20    (FPT) 
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I4.O.  transforms  averaged  per  trace,  transform  length  20i+8 
Displayed  bandv;idth  l8,75  hz,  0,125  hz  resolution 

Figure  2-20  (Mod-FF"^"* 
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14.0  transforms  averaged  per  trace,  transform  length  201^8 
Displayed  bandwidth  18.75  hz,  0.12^  hz  resolution 


S-fM-tBU.^ST     TFT 

S/N'"~30  DB         Figure  2-21  (FFT) 
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I4.O  transforms  averaged  per  trace,  transform  length  20i|.8 
Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 


Figure  2-21  (Mod-FFT) 
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FFT 
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Modified  FFT 

10  hz  signal,  S/N  =  -12  db  +  burst  noise 

40  transforms  averaged  per  trace,  transform  length  204  8 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-22 
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Modified  FFT 

10  hz  signal,  S/N  =  -18  db  +  burst  noise 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-23 
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FFT 


Modified  FFT 

10  hz  signal,  S/N  =  -24  db  +  burst  noise 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-24 
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FFT 


Modified  FFT 


10  hz  signal,  S/N  =  -24  db  +  burst  noise 

80  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-25 


90 


FFT 


Modified  FFT 

10  hz  signal,  S/N  =  -24  db  +  burst  noise 
80  transforms  averaged  per  trace,  transform  length  204  8 
Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 
Spectrum  squared 

Figure  2-26 
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Modified  FFT 

Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Figure  2-27 


92 


FFT 


Modified  FFT 


Displayed  bandwidth  18.75  hz,  0.125  hz  resolution 

Spectrum  squared 


Figure  2-28 
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S/N  out  (db) 


S/N  in 
(db) 

FFT 

Mod-FFT 

FFTxFFT 

FFTx 
Mod-FF' 

-12 

13.770 

22.354 

16.342 

23.711 

-18 

7.729 

14.389 

8.671 

15.133 

-24       8.747      8.360     6.450     8.863 
-30        5.109      5.515      3.551      5.910 

Table  IV 

Average  S/N  out  (db)  vs.  S/N  in  (db)  for  the  Average  of  40 

Transforms  for  the  case  of  the  10  hz  signal  plus  Gaussian 

noise  plus  burst  noise.   S/N  out  computed  as: 

(F  -y  )2 
^  s   n 


2a  ^ 
n 

where  F   is  the  magnitude  of  the  signal  coefficient  and  y 
s  '^  '^  n 

and  a   are  computed  from  25  noise  coefficients  on  either 

side  of  F  . 
s 
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S/N  out  (db) 

^/!Jun"       FFT      Mod-FFT    FFTxFFT    ./^'^^t^m 
(db)  Mod-FFT 

-12  16.154  26.606  18.731  26.945 

-18  11.433  16.341  12.855  17.121 

-24  9.332  10.979  6.713  10.061 

-30  6.875  8.094  5.130  7.961 

Table  V 

Average  S/N  out  (db)  vs.  S/N  in  (db)  for  the  Average  of  80 

Transforms  for  the  case  of  the  10  hz  signal  plus  Gaussian 

noise  plus  burst  noise.   S/N  out  computed  as: 

(F  -M    )2 
^  s   n 


where  F   is  the  magnitude  of  the  signal  coefficient  and  y 
s  ^  n 

and  a   are  computed  from  25  noise  coefficients  on  either 

side  of  F  . 
s 
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E.   UNDERWATER  ACOUSTIC  DATA 

The  final  data  to  be  analysed  was  underwater  acoustic 
data.   This  data  provides  a  great  variety  of  signals  from 
different  sources  with  different  strengths  all  transmitted 
simultaneously.   It  is  a  challenging  test  of  the  signal  proc- 
essor's ability  to  resolve  the  many  different  signals  and  to 
perceive  weak  signals  in  the  presence  of  strong  ones.   Only 
selected  portions  of  the  data  have  been  presented  in  this 
section. 

The  data  was  filtered  at  200  hz  and  sampled  at  512  sam- 
ples per  second.   Again,  2048  point  transforms  were  employed 
resulting  in  a  frequency  resolution  of  J  hz. 

Each  trace  in  the  pictures  is  the  result  of  averaging 
40  transforms.   The  FFT  spectra  are  in  the  top  photographs 
in  each  figure  and  the  modified  FFT  spectra  are  in  the  lower 
photographs.   The  time  frame  and  frequency  bandwidth  are  the 
same  for  each  photograph  in  a  given  figure.   The  total  dis- 
played frequency  bandwidth  in  each  trace  is  37.5  hz. 

Figures  2-29  and  2-30  each  depict  the  same  frequency 
bandwidth  but  at  different  times.   The  main  distinction  is 
the  modified  FFT's  ability  to  detect  the  highlighted  signal 
in  Figure  2-29  before  the  FFT  and  to  make  it  clearly  recog- 
nizable.  The  large  signals  to  the  left  and  right  of  the 
highlighted  signal  are  enhanced  better  by  the  FFT.   There 
are  also  significant  differences  in  the  general  appearance 
of  the  two  spectra.   Figure  2-30  is  presented  to  show  that 
at  a  later  time  the  FFT  presented  the  same  signal  as  well  as 
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or  better  than  the  modified  FFT.   The  signal  to  the  left  of 
the  highlighted  signal  seems  to  be  presented  equally  well  by 
both  spectra. 

Figures  2-31  and  2-32  show  a  relatively  weak  signal. 
Close  comparison  of  the  traces  in  the  spectra  indicate  that 
in  some  cases  the  FFT  outperforms  the  modified  FFT  and  in 
others  the  reverse  is  true.   The  bottom  trace  of  Figure  2-32 
provides  a  dramatic  case  in  which  the  modified  FFT  produces 
the  signal  better  than  the  FFT. 

Figures  2-33  through  2-35  provide  an  example  in  which 
the  two  spectra  behave  almost  identically  except  for  the 
last  trace  in  Figure  2-35  in  which  an  effect  like  the  one 
described  for  Figure  2-32  above  occurs. 

Figure  2-36  again  depicts  a  signal  that  has  been  en- 
hanced better  by  the  modified  FFT  than  by  the  FFT.   Also  two 
distinct  signals  are  present  to  the  right  of  the  highlighted 
signal  in  the  modified  FFT  whereas  only  one  is  present  in 
the  FFT. 

Figures  2-37  through  2-39  all  show  the  same  signal  at 
different  times.   In  this  case  the  two  spectra  are  again 
quite  similar.   The  main  reason  for  presenting  these  spectra 
is  the  modified  FFT's  ability  to  retain  the  signal  in  Figure 
2-39  after  it  has  been  lost  by  the  FFT. 
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FFT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-29 
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FFT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-30 
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FFT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-31 
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FFT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-32 
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FFT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  204  8 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-33 
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PPT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-34 
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FFT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-35 
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FFT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-36 
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FFT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  204  8 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-37 
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FFT 


Modified  FFT 

Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-38 
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FFT 


Modified  FFT 


Underwater  Accoustic  Data 

40  transforms  averaged  per  trace,  transform  length  2048 

Displayed  bandwidth  37.5  hz,  0.25  hz  resolution 

Figure  2-39 
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IV.   PART  3  SUMMARY  OF  EXPERIMENTAL  RESULTS 

The  data  analysed  in  this  paper  were  basically  of  three 
types.   A  signal  or  signals  with  additive  Gaussian  noise,  a 
signal  with  additive  Gaussian  noise  plus  bursts  of  Gaussian 
noise  spaced  at  random  intervals  with  random  duration,  and 
finally  some  underwater  acoustic  data. 

In  Section  III-B  of  this  investigation,  a  single  sinu- 
soid with  additive  Gaussian  noise  was  processed  with  the  FFT 
and  the  modified  FFT.   Figures  2-1  through  2-9  indicated 
that  the  resulting  spectrums  from  the  two  transformations 
are  quite  similar.   The  FFT  seems  to  have  smoothed  the  noise 
slightly  better  than  the  modified  FFT.   In  general  the  FFT 
seems  to  provide  the  better  processing  gain,  but  upon  close 
scrutiny  of  the  spectra  it  is  possible  to  find  some  traces 
in  which  the  modified  FFT  seems  to  produce  a  greater  signal 
to  noise  ratio.   The  computed  S/N  output  ratios  in  Table  I 
seem  to  agree  quite  well  with  what  is  observed  in  the  photo- 
graphs.  The  modified  FFT  seems  to  outperform  the  FFT  at  the 
large  S/N  ratios  and  the  FFT  dominates  the  modified  FFT  at 
the  small  signal  to  noise  ratios. 

Section  III-C  of  this  paper  dealt  with  multiple  signal 
data  with  additive  Gaussian  noise.   Figures  2-10  through 
2-14  depict  the  output  spectra  and  again  it  is  noted  that 
the  spectra  are  quite  similar  in  overall  appearance,  with 
neither  transformation  appearing  to  have  a  distinct  advan- 
tage over  the  other.   The  most  interesting  results  are  in 
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Figures  2-13,  where  in  the  FFT  spectra  the  -16  db  signal  is 
consistently  smaller  than  the  -14  db  signal.   In  the  modi- 
fied FFT  spectra  the  -16  db  signal  is  nearly  as  large  as  the 
-14  db  signal.   From  the  -12  db  signal  on  down  it  is  doubt- 
ful that  one  could  recognize  the  signal  presented  in  the 
figures.   Longer  transforms  or  the  averaging  of  more  trans- 
forms would  be  required  to  enhance  the  signal  to  noise  ra- 
tios.  Averaging  twice  as  many  transforms  theoretically 
increases  the  S/N  output  by  3  db  and  this  is  verified  in 
Tables  II  and  III. 

These  tables  tend  to  indicate  that  the  FFT  is  the  better 
transformation  at  high  S/N  ratios  and  the  modified  FFT  is 
better  at  low  S/N  ratios.   For  reasons  pointed  out  in  Section 
III-C  the  -18  db  signal  should  not  be  used  in  the  analysis. 

Some  of  the  apparent  inconsistencies  in  these  tables  are 
explained  by  the  fact  that  different  noise  coefficients  were 
used  in  computing  the  noise  statistics  for  each  signal  as 
mentioned  in  Section  III-C.  There  is  no  known  explanation 
for  the  inconsistency  of  the  process  gain  for  the  -6  db  sig- 
nal. A  search  through  the  adjacent  coefficients  showed  that 
this  coefficient  provided  the  best  signal  to  noise  output  in 
the  vicinity  of  the  signal  for  the  FFT. 

Another  interesting  aspect  of  Tables  II  and  III  is  that 
in  going  from  a  40  transform  average  to  an  80  transform  aver- 
age the  modified  FFT  S/N  ratios  are  improved  more  than  the 
FFT  S/N  ratios. 
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Section  III-D  presents  data  in  which  bursts  of  noise 
were  added  to  the  signal  in  addition  to  continuous  Gaussian 
noise.   In  this  case  the  modified  FFT  consistently  outper- 
formed the  FFT.   The  one  exception  was  in  the  case  of  the 
-24  db  signal  where  40  transforms  were  averaged.   When  80 
transforms  were  averaged  the  modified  FFT  outperformed  the 
FFT  in  every  case.   Figures  2-25  through  2-28  show  the  re- 
sults of  the  80  transform  average  in  which  it  is  much  easier 
to  believe  in  the  existence  of  a  signal  in  the  case  of  the 
modified  FFT  than  in  the  FFT. 

Tables  IV  and  V  present  the  S/N  ratios  for  the  spectra. 
Little  needs  to  be  said  in  this  case  since  the  results  seem 
to  be  consistent  with  the  observations  made  in  Figures  2-18 
through  2-28. 

Section  III-E  of  this  paper  presents  several  signals  en- 
countered in  the  transformation  of  some  actual  underwater 
acoustic  data,  an  area  in  which  signal  processing  finds  a 
very  practical  use.   The  manner  of  interpretation  of  the  data 
is  very  subjective,  as  is  the  manner  in  which  it  is  normally 
analyzed.   In  the  signals  that  were  presented  the  modified 
FFT  usually  did  as  well  or  better  than  the  FFT  in  making  sig- 
nal presence  easily  detectable.   During  all  of  the  research 
that  was  done  it  was  found  that  both  transformations  gener- 
ally provided  similar  results  for  this  type  of  data.   Ex- 
tremely strong  signals  were  generally  enhanced  more  by  the 
FFT.   Medium  to  weak  signals  appeared  to  be  about  the  same 
in  overall  detectability  for  both  transl'ormation.s .   Looking 
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specifically  from  trace  to  trace  differences  were  definitely 
noted  and  no  one  particular  transformation  proved  to  be  con- 
sistently better.   On  some  occasions  the  modified  FFT  pro- 
duced signals  that  the  FFT  could  not,  even  though  later  on 
the  signal  was  produced  by  the  FFT.   In  the  case  of  Figure 
2-29  an  analyst  confirmed  the  high  probability  of  a  signal 
being  present  at  the  spotlighted  frequency.   In  Figure  2-39 
the  ability  of  the  modified  FFT  to  retain  the  signal  longer 
than  the  FFT  is  another  strong  point  in  favor  of  the  modi- 
fied FFT  but  an  unusual  occurrence.   The  nonlinearity  of  the 
modified  FFT  with  the  transformation  being  a  function  of  the 
input  data  itself  has  made  it  difficult  to  draw  generalities 
about  its  performance  except  possibly  in  the  case  of  data 
similar  to  that  presented  in  Section  III-D  of  this  paper. 
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V.   PART  4  CONCLUSIONS 

The  data  analyzed  in  this  paper  is  far  from  conclusive 
or  absolute  in  terms  of  allowing  any  hard  and  fast  conclu- 
sions to  be  drawn.   The  process  does  seem  to  have  some  mer- 
its, especially  in  the  case  of  the  signal  plus  noise  plus 
burst  noise.   Whereas  this  may  not  appear  to  be  a  very  prac- 
tical type  of  data,  it  is  felt  that  somewhere  between  this 
case  and  the  case  of  stationary  Gaussian  noise  a  practical 
situation  exists  in  which  this  transformation  may  be  very 
useful.   Somewhere  where  strong  transients  exist,  whether 
they  are  switching  noises,  or  the  explosion  of  underwater 
sonic  sources,  this  transformation  may  be  useful  if  an  under- 
lying continuous  signal  is  of  interest. 

Unexplained  are  the  instances  in  Section  III-E  of  this 
paper  where  the  modified  FFT  was  able  to  bring  out  signals 
significantly  better  than  the  FFT,  or  in  one  case  to  retain 
the  signal  for  a  longer  period  of  time  than  the  FFT. 

On  the  other  hand,  the  FFT  seems  better  able  to  distin- 
guish between  signal  and  noise  when  the  noise  is  stationary 
Gaussian. 

The  good  characteristics  of  the  transformation  do  not 
seem  to  be  dependent  on  the  signal  frequency.   In  Section 
III-E  the  signals  that  were  much  more  evident  in  the  modi- 
fied FFT  than  in  the  FFT,  appeared  at  both  ends  of  the 
spectrum  as  well  as  in  the  middle  of  the  spectrum.   This  ef- 
fect was  not  limited  to  one  data  tape  but  appeared  on  two 
different  tapes  from  different  sources. 
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The  main  drawback  to  implementing  the  modified  FFT  as  it 
is  presently  being  utilized  is  the  length  of  time  required 
to  perform  the  transformations,  approximately  three  to  four 
times  as  long  as  a  similar  FFT  transformation  for  1024  points. 
However  if  parallel  processing  were  available  the  modified 
FFT  could  be  performed  twice  as  fast  as  the  FFT  since  each 
pass  only  operates  on  half  of  the  coefficients. 

In  general  there  does  not  seem  to  be  enough  difference 
between  the  FFT  and  modified  FFT  spectra  to  warrant  its  use 
in  preference  to  the  FFT.   If  it  is  to  be  used  it  would 
seem  appropriate  to  use  it  when  a  signal  is  suspected  and 
the  FFT  does  not  reveal  it.   Also,  if  the  data  were  known  to 
have  characteristics  similar  to  that  of  the  signal  plus  noise 
plus  burst  noise  it  would  warrant  strong  consideration. 
One  other  asset  of  the  transformation  that  may  prove  impor- 
tant is  the  inability  of  any  given  output  coefficient  to 
exceed  some  maximum  value,  determined  by  the  transform  length, 

With  so  many  variables,  such  as  transform  length,  fre- 
quency resolution,  windowing,  signal  stability,  and  differ- 
ing types  of  background  noise  and  signals,  much  more 
investigation  will  be  required  to  determine  the  optimum 
utilization  of  this  transformation.   A  mathematical  and 
statistical  analysis  of  the  modified  FFT  may  give  further 
insight  into  its  best  application. 
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APPENDIX  A 
FAST  FOURIER  TRANSFORM  " 

The  Fast  Fourier  Transform  is  an  algorithm  for  computing 
the  "Discrete  Fourier  Transform"  which  evaluates  the  Fourier 
Coefficients  of  the  Fourier  Series  for  a  discretely  sampled 
n-wwaveform. 

The  term  Fast  Fourier  Transform  is  somewhat  misleading 
because  the  FFT  is  applied  to  the  transformation  of  periodic 
functions  which  do  not  necessarily  satisfy  the  condition: 

00 

/  |f(t)|  dt  <  - 

— oo 

which  is  required  for  f(t)  (in  general)  if  it  is  to  have  a 
Fourier  Transform. 

Any  discrete  analysis  of  a  function  f(t)  must  of  neces- 
sity be  finite  in  length,  say  from  t=0  to  t=T.   What  happens 
before  t=0  and  after  t=T  is  in  general  not  known,  or  at 
least  the  time  record  in  question  gives  no  indication  of 
what  was  occurring  prior  to  t=0  or  what  is  going  to  occur 
after  t=T,  unless  the  function  is  periodic  and  the  time  re- 
cord contains  at  least  one  whole  period  of  the  waveform. 
AVhen  dealing  with  the  FFT  the  implied  assumption  is  that  the 
record  to  be  analyzed  is  one  period  of  a  periodic  function 
and  that  the  information  contained  in  this  period  completely 
describes  the  function  from  t^-oo  to  t=+o°. 

Periodic  functions  are  more  appropriately  described  by 
their  Fourier  Series  expansion,  rather  than  by  the  Fourier 
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Transform.   The  term  Fast  Fourier  Transform  may  be  more  de- 
scriptive if  it  were  renamed  the  Fast  Fourier  Series. 

A  waveform  f(t),  periodic  in  t  with  period  T  satisfying 
f(t)=f(t+T)  for  all  t  with  a  finite  number  of  discontinuities 
on  the  interval  (0<_t_<T)  can  be  expanded  in  a  Fourier  Series 
as: 

oo 

f(t)  =  I      F(n)exp[(  j2TTnt)/T] 
n=-oo 

where  F(n),  the  Fourier  Coefficients  are  defined  as: 

1  ^ 
F(n)  =  i  /  f(t)exp[(-j2TTnt)/T]dt 

^  0 

The  above  representation  of  the  Fourier  Coefficients  is 
for  the  continuous  case.   In  order  to  determine  the  coeffi- 
cients for  the  discrete  case  the  integral  must  be  approxi- 
mated by  a  summation.   Assuming  the  waveform  has  been  sampled 
N  times  in  T  seconds  and  that  the  samples  are  spaced  AT  sec- 
onds apart,  then, 

T=NAT 
and 

t=kAT   (k=0,l,2, . . . ,N-1) 
then  the  sampled  version  of  the  Fourier  Coefficients  can  be 

written : 

N-1 
F(n)    =  ^    S      f(kAT)exp[(-j27TkAT)/T]AT 

N-1 
=i      Z      f(kAT)exp[(-j2TTnkT)/(NT)] 
^^    k=0 

N-1 
=^      E      f(kAT)exp[(-j2Trnk)/N]       (  n=0  , 1 ,  2  ,  .  .  .  ,N-1 ) 
^    k-0 
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The  factor  1/N  is  a  scaling  factor,  useful  for  obtaining 
the  original  amplitude  of  the  input  waveform  when  f(t)  is 
expressed  in  its  Fourier  Series.   It  has  no  effect  on  the 
relative  amplitudes  of  the  Fourier  Coefficients. 

F(n)  is  the  amount  of  frequency  n/T  contained  in  f(t). 

Abbreviating  f(kAT)  as  f(k),  and  omitting  the  scaling  factor 

1/N  for  convenience,  then, 

N-1 
F(n)  =   E   f(k)exp[(-j27Tnk)/N] 
k=0 

N-1 
=      Z      f(k){cos[(2TTnk)/N]-jsin[(2TTnk)/N]} 
k=0 

N-1  N-1 

=      E      f  (k)cos[(2unk)/N]-j       E      f  (k)sin  [  (2TTnk) /N] 
k=0  k=0 

which  is  simply  a  discrete  correlation  of  f(t)  with  N  sinu- 
soids at  the  frequencies  n/T  since, 

2-n-nk  ^  2TTtn  ^  2iTntAT  ^  2TTnt     ,     H  =   Af 
N      NAT      TAT       T      wnere  ^    nAi 

F(0)  is  then  the  amount  of  D.  C.  component  contained  in  f(t) 

F(l),  the  amount  of  frequency  1/T  =  Af  contained  in  f(t), 

F(2),  the  amount  of  frequency  2/T  =  2Af  contained  in  f(t), 

etc.,  up  to  n=N-l  for  a  complex  waveform  f(t). 

If  f(t)  is  real  valued  the  Fourier  Coefficient  F(N/2-m) 

=  F  (N/2+m)  (where  *  denotes  the  complex  conjugate). 

Proof : 

1   N-^ 
F(n)    =   M      E      f(k){cos[(2Tink)/N]-jsin[(2Tink)/N]  } 

^    k=0 
let   F(n)    =   R(n)-jl(n),    then. 
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1 N-1 

R(n)    =   Re[F(n)]    =  ^      T.      f  (k)cos  [  (2Trnk)/N] 

^   k=0 

for   n=N/2,    R(N/2)=  ^      Z      f (k)cos(TTk)=  ^      Z    (-1)    f(k) 

^   k=0  ^   k=0 

now:    cos[27rk(N/2-m)/N]    =    cos  [  (N7rk-2TTkm)/N] 

=   cos[(NTTk)/N]cos[(2Trkm)/N]    +    sin  [  (NTTk)/N]  sin  [  (2TTkm)/N] 
where   the   second   term   is   zero   for   integer  k. 

=   cos[(NkTT   +   2iTkni)/N]    =   cos  [2TTk(N/2   +  m)/N] 
therefore, 

R(N/2-ra)  =  R(n/2+m)   (m=0, 1, 2 , . . . ,N/2-l) 
or  the  real  part  of  F(n)  is  symmetric  about  n=N/2. 

A  similar  proof  shows  that  the  imaginary  part  of  F(n)  is 
antisymmetric  about  n=N/2,  i.e., 

I(N/2-m)  =  -I(N/2+m)   (m=0, 1 , 2 , . . . , N/2-1) 


■nil 


Figure  Al-1 


That  is  F(N/2-m)  =  F  (N/2+m),  implying  that  all  of  the 
spectral  information  is  contained  in  the  first  N/2  coeffi- 
cients, the  second  half  of  the  spectrum  containing  redundant 
information.   This  is  intuitively  obvious  if  one  considers 
that  the  cosine  is  an  even  function  and  the  sine  is  an  odd 

function. 
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To  more  readily  visualize  the  correspondence  between  the 
Fourier  Transform  of  a  continuous  and  a  sampled  waveform, 
consider  a  continuous  waveform  f(t),  and  its  samples  counter- 
part f  (t)  (where  the  k  denotes  the  sampled  version  of  f(t)). 
Then  f|^(t)  can  be  represented  as: 

fj^(t)  =  f(k)  =   E   f(t)6(t-kAT) 

k=-<» 

CO 

let  d(t)  =   E   6(t-kAT) 
k=-oo 

D(w)  ^d(t) 
F(w)  <^  f  (t) 
Fj^(w)  ^fj^(t) 

where  <^  implies  Fourier  Transform  pairs. 
Then  F  (w)  =  D(w)*F(w) 

Since  {6(t-kAT)}  is  a  train  of  impulse  functions  periodic 
in  t  with  period  AT,  the  Fourier  Series  representation  is 
used  to  describe  this  function,  that  is: 

oo 

d(t)    =      E      D(n)exp(jwnt)      where  w=2TTf   =    2tt/AT 
n=-oo 

and 

^    AT 
D(n)    =  ^   f   d(t)exp(-jwnt)dt 

then, 

^     AT   CO 
D(n)    =  ~   f      E      6(t-kAT)exp(-jwnt)dt 

k=-oo 


AT    0    ^=. 


^       CO         AT 

-~      L         f    6(t-kAT)exp(-jwnt)dt 

A  i     1  _  0 

k=-oo 
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1  1 

D(n)    =  -^   ^      exp(-jwntkAT)    =  t^   2      exp  [(- j2TrnkAT)/AT  ] 
k=-oo  k=-oo 

OO 

1  1 

~    AT    ^      exp(-j27Tnk)    -  — =      for    all    integer   n,    k 
k=-oo 


Now, 


d(t)    =    6(t-kAT)    =  -r|   Z      exp(jwnt) 

n=-oo 

then: 

CO  -1         °° 

Fj^(w)    =      •'"[^(t^AT    ^      exp(jwnt)]exp(-jwt)dt 

—  OO  j;^  =  _oo 

OO  OO 

^    AT    ^  /[f(t)exp    ( jwnt) ]exp(-jwt)dt 

ri=— OO  — OO 

OO  CO 

^    AT    ^  /   f (t)exp[j(wn-w)t ]dt 

J^  =  — CO       — OO 
CO  OO 

^    AT    ^         /   f  (t)exp[-3(wn-w)t  ]dt 

j;i=— OO     —CO 
1     ~ 

^    AT    ^      F(w-wn) 
n=-oo 

The  discrete  Fourier  Transform: 

F(n)  =^  ^  f  (kAT)exp(-j2TTnk/N)  (n-0, 1 ,  2  ,  .  .  .  ,N-1) 

can  be  visualized  in  various  ways.   As  previously  mentioned, 
it  can  be  considered  a  discrete  correlation  of  the  sampled 
waveform  f(kAT)  with  a  series  of  sinusoids  at  the  frequencies 
n/T  (n=0,l,2, . . . ,N-1)  where  T=NAT  is  the  length  of  the  input 
waveform  to  be  transformed,  N  is  the  total  number  of  samples, 
and  AT  is  the  time  between  samples. 
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The  time  T  then  determines  the  frequency  resolution  of 
the  transformation,  that  is,  l/T=Af  or  the  difference  in 
frequency  from  one  correlating  sinusoid  to  the  next.   N  and 
T  determine  the  highest  frequency  of  the  correlating  sinu- 
soids.  The  ratio  of  N/T  is  the  sampling  rate,  i.e.,  N  sam- 
ples in  T  seconds.   For  real  valued  input  functions,  the 
highest  correlating  frequency  is  then  (N/2-l)/T  or  half  of 
the  sampling  rate.   Figure  Al-2  shows  that  the  waveform  must 
be  sampled  at  a  rate  at  least  twice  as  high  as  the  highest 
frequency  present  in  f(t)  or  the  spectra  of  F(w)  will  over- 
lap.  This  is  called  aliasing. 

A  second  visualization  is  simply  that  of  a  matrix  opera- 
tion.  Expressing  the  discrete  Fourier  Transform  in  matrix 
form,  omitting  the  scaling  term  1/N,  and  defining  V/  as  exp 
(-J21T/N)  then: 

[F(n)]  =  [W''^][f(k)] 
An  expansion  of  these  matrices  for  eight  points  is  given  in 
Figure  Al-3a. 

Figure  Al-3b  shows  the  matrix  rewritten  with  the  follow- 
ing substitutions: 

^^nk  ^  ^^,nkModN   g^j^^^  ^^^  is  periodic  with  period  N. 
Also,    W"  =  exp(O)  =  1  =  -W" 
Wi  =  exp(-J7r/4)  =  -W^ 
W2  =  exp(-JTr/2)  =  -W^ 
W    =   exp(-j37T/4)  =  -W^ 
The  matrix  in  Figure  Al-3b  can  be  expanded  into  three 
matrices  as  shown  in  Figure  Al-4. 
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Discrete  Fourier  Transform  in  matrix  form 
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1   1 
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1  " 

f(o) 

F(1) 

1    W^ 

w2 
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1  -w^ 
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f(1) 

F(2) 
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-1 

-w2  - 

1    W2 

-1 
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f(2) 
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W2 
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Matrix  in  (a)  rewritten  usin  the  substitution 

,,nk  _  ,,nkModN 
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Figure  A1-3 
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It  is  this  matrix  expansion,  along  with  the  periodic 

nk 
properties  of  W    (the  correlating  sinusoids),  that  gives 

rise  to  the  Fast  Fourier  Transform.   In  general,  if  a  sam- 

M 
pled  input  waveform  has  N  samples  where  N=2   and  M  is  a  pos- 
itive integer,  then  the  NxN  matrix  of  operators  (the  matrix 
of  W's)  can  be  expanded  into  a  product  of  M  matrices  with 
the  stipulation  that  the  input  vector  is  bit  reversed  (note 
the  reordering  of  the  input  vector).   This  stipulation  is 
not  a  requirement  but  it  is  the  method  that  is  used  through- 
out this  study.   The  process  of  bit  reversal  allows  the  ma- 
trix multiplications  to  be  performed  one  after  the  other 
without  a  requirement  to  save  any  intermediate  results.   The 
matrix  expansion  can  be  performed  in  other  ways  such  that 
the  input  vector  is  left  in  its  natural  order  and  the  output 
vector  is  in  bit  reversed  order,  or  possibly  the  odd  harmon- 
ics of  Af  are  grouped  together  and  the  even  harmonics  are 
grouped  together.   It  is  also  possible  to  perform  the  opera- 
tions without  bit  reversing  the  input  or  the  output  vectors 
but  then  a  requirement  exists  for  the  retention  of  inter- 
mediate results. 

Bit  reversal  is  simply  the  reversing  of  the  order  of  the 
digits  that  make  up  a  number  in  binary  arithmetic.   In  the 
above  situation  N-8,  and  M=3.   For  example,  110  (6)  becomes 
Oil  (3)  when  its  bits  are  reversed.   The  effect  of  bit  re- 
versal is  to  place  each  element  in  the  input  vector  next  to 
an  element  that  is  N/2  elements  away  from  it,  i.e.,  N/2-4  in 
this  case  so  element  1  is  placed  next  to  element  5,  2  next 

to  6,  etc. 
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To  go  one  step  further,  the  matrices  of  Figure  Al-4  can 
be  further  reduced  by  using  a  tree  graph  or  what  is  commonly 
called  an  "FFT  butterfly".   The  expanded  matrices  have  many 
zero  elements,  in  fact  each  matrix  contains  only  2N  elements 
that  are  non-zero  and  a  large  portion  of  these  are  1 ' s  or 
-I's.   The  original  Discrete  Fourier  Transform  required  N^ 
complex  operations  in  its  execution;  the  expanded  matrices 
of  the  Fast  Fourier  Transform  only  have  Nlog2N  operations 
and  thus  the  term  "Fast". 

The  matrices  have  been  rewritten  in  butterfly  format  in 
Figure  Al-5.   The  column  labeled  "0"  is  the  input  vector. 
The  column  labeled  "1"  is  the  bit  reversed  form  of  the  input 
vector.   The  column  labeled  "2"  is  the  results  of  the  first 
matrix  multiplication,  or  the  results  after  the  first  pass 
through  the  algorithm.   The  number  2  is  used  to  describe 
this  column  since  two  of  the  original  elements  make  up  each 
of  the  elements  at  this  point  in  the  algorithm.   The  column 
labeled  "4"  is  the  result  of  the  second  matrix  multiplica- 
tion, or  the  result  after  the  second  pass  through  the  algo- 
rithm.  As  above  the  number  4  is  used  to  describe  this 
column  since  four  of  the  original  elements  make  up  each  of 
the  elements  at  this  point  in  the  algorithm.   The  nomencla- 
ture is  similar  for  each  of  the  columns  in  the  butterfly. 

When  two  arrows  meet  at  a  circle  with  a  +  or  -  in  it, 
then  addition  or  subtraction  is  to  take  place.   For  example, 
in  the  column  labeled  2  element  fo  is  made  up  of  elements 
fo  +  fi  from  column  1;  element  fj  in  column  2  is  made  up  of 
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f-,  =  F(1) 


b)f2  =  F(2) 


^3  =  F(3) 
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FPT  butterfly  or  flow  diagram  for  an  8  point 
transformation 


Figure  A1-5 
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elements  fo  -  fi  from  column  1.   The  circles  that  contain  an 
argument  like  1,  2,  or  3  imply  a  multiplication  of  the  pre- 
ceding element  by  W\  W^ ,  or  W^  before  an  addition  or  sub- 
traction is  performed.   For  example,  in  column  4  element 
fo  is  made  up  of  elements  fo  +  f2  from  coliimn  2,  element 
fi  is  made  up  of  elements  fi  +  fsW^  from  column  2,  element 
iz    is  made  up  of  elements  fo  -  fa  from  column  2,  element 
fs  is  made  up  of  elements  fi  -  fsW^  from  column  2. 

The  butterfly  is  identical  to  the  matrix  approach  in 
that  each  pass  through  the  butterfly  is  identical  to  perform- 
ing one  of  the  matrix  multiplications  in  the  matrix  repre- 
sentation of  the  FFT.   The  butterfly  condenses  the  matrix 
approach  by  elimination  of  the  zeros  in  each  matrix  leaving 
only  the  necessary  operations. 
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